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Abstract. Peak stratospheric chlorofluorocarbon (CFC) and 
other ozone depleting substance (ODS) concentrations were 
reached in the mid- to late 1990s. Detection and attribution 
of the expected recovery of the stratospheric ozone layer in 
an atmosphere with reduced ODSs as well as efforts to un- 
derstand the evolution of stratospheric ozone in the presence 
of increasing greenhouse gases are key current research top- 
ics. These require a critical examination of the ozone changes 
with an accurate knowledge of the spatial (geographical and 
vertical) and temporal ozone response. For such an exami- 
nation, it is vital that the quality of the measurements used 
be as high as possible and measurement uncertainties well 
quantified. 

In preparation for the 2014 United Nations Environ- 
ment Programme (UNEP)/World Meteorological Organiza- 
tion (WMO) Scientific Assessment of Ozone Depletion, the 
SPARC/IO 3 C/IGACO-O 3 /NDACC (SI 1 2 N) Initiative was de- 
signed to study and document changes in the global ozone 
profile distribution. This requires assessing long-term ozone 
profile data sets in regards to measurement stability and 
uncertainty characteristics. The ultimate goal is to estab- 
lish suitability for estimating long-term ozone trends to con- 
tribute to ozone recovery studies. Some of the data sets have 
been improved as part of this initiative with updated versions 
now available. 

This summary presents an overview of stratospheric ozone 
profile measurement data sets (ground and satellite based) 
available for ozone recovery studies. Here we document 
measurement techniques, spatial and temporal coverage, ver- 
tical resolution, native units and measurement uncertainties. 
In addition, the latest data versions are briefly described (in- 
cluding data version updates as well as detailing multiple 
retrievals when available for a given satellite instrument). 
Archive location information for each data set is also given. 

1 Introduction 

As man-made ozone depleting substances (ODSs) decline in 
the stratosphere in response to the 1987 Montreal Protocol 
and its subsequent amendments and adjustments, the ozone 
layer is expected to recover globally. A small increase in 
stratospheric ozone abundance over northern mid-latitudes 
was first reported in the 2010 WMO ozone assessment 


(WMO, 2011). Although supported by a variety of ground- 
based, balloon-born and satellite observations, statistical 
analysis could not attribute this change exclusively to the de- 
cline of ODSs due to observational uncertainty, geophysical 
and dynamical variability and changes in ozone due to in- 
creasing greenhouse gases that cool the stratosphere. Global 
total ozone remained significantly lower than during the early 
1980s, and surface ultraviolet radiation measurements were 
too variable to strongly support ozone recovery. Indications 
of the onset of stratospheric ozone recovery have been noted 
since the publication of the last ozone assessment. There are 
several recent publications (e.g. Maeder et al., 2010; Nair 
et al., 2013; Kuttippurath et al., 2013) supporting the state- 
ment that the second stage of ozone recovery is ongoing 
(as defined by the WMO, 2007, statement that “occurrence 
of statistically significant increases in ozone above previous 
minimum values due to declining equivalent effective strato- 
spheric chlorine”). Other publications show that ozone re- 
covery in Antarctica should be detectable within the next 
decade (Newman et al., 2006; Hassler et al., 2011). The ad- 
ditional four years of data available since analysis was done 
for the 2010 WMO ozone assessment as well as improve- 
ments in the homogenization of ozone data records should 
improve the signal-to-noise ratio in the ozone time series 
making trend analyses results more robust. 

“SPARC Report No.l - Trends in the vertical distribution 
of ozone” (SPARC, 1998) provided a comprehensive sum- 
mary of available data and methods for trend analysis at that 
time. Its main objectives were to understand the limitations 
of the data, to assess the consistency between various data 
records and create a reference database for the analysis of 
long-term ozone changes as a function of altitude. One of the 
main reasons for this approach was to avoid misinterpreting 
instrument drifts as actual ozone changes. 

Trend analysis is made more complex by the addition of 
new satellite records (e.g. Aura, Envisat, Odin, and SCISAT), 
where each covers a shorter time period than the tradition- 
ally used long-term records (e.g. ozonesondes, Umkehr, li- 
dars, SAGE, HALOE). Adding to the problem, some of 
the long-term platforms have ceased taking measurements 
(e.g. SAGE, HALOE, instruments on Envisat), while records 
from other satellite platforms are difficult to combine due to 
drifting orbits (SBUV/2 on NOAA satellites) and different 
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instrument generations (e.g. moving from SBUV fixed wave- 
length retrievals to Aura OMI hyper-spectral retrievals). For 
a full long-term analysis of stratospheric ozone throughout 
the decline due to increases in ODSs and subsequent re- 
covery, multiple instrument records need to be considered. 
Thus, several different data sources have to be merged to al- 
low a complete view of past ozone changes. Recent work 
has shown that ozone depletion and greenhouse gas (GHG) 
increases impact dynamical quantities in an additive man- 
ner, with ozone recovery possibly cancelling atmospheric 
circulation effects associated with increasing GHGs. (e.g. 
IPCC/TEAP, 2005; McLandress et al., 2010; Polvani et al., 
2011). Clearly, an understanding of atmospheric transport, 
circulation and wave breaking patterns is important for our 
understanding of ozone recovery. 

The consideration of uncertainties and artifacts is crucial 
in ozone trend analyses, particularly when the ozone recov- 
ery signal is small as compared to natural geophysical vari- 
ability. As methods for trend and attribution analysis become 
more advanced, the uncertainties and instrumental artifacts 
buried in currently available data sets start to play a bigger 
role. Therefore, detailed information about measurement un- 
certainties, data jumps due to changes in instrumentation and 
drifts during long-term operation of a given instrument are 
essential. 

There have been several attempts to produce histori- 
cal global ozone profile data sets; these include satellite- 
ozonesonde combinations starting in 1979 when the first 
global satellite measurements became available (e.g. Randel 
and Wu, 2007; Bodeker et al., 2013), global ozone data sets 
that extend values to 2100 using model and fitting techniques 
(Hassler et al., 2009; Cionni et al., 2011) for use in global 
climate models without stratospheric ozone chemistry and 
also data sets that extend the historical record backward to 
1900 (Bronnimann et al., 2013). Issues of biases between in- 
struments, sampling, noise, gaps and differences in vertical 
registration (pressure or altitude) have to be addressed when 
creating combined data records to ensure no artificial trends 
are introduced. 

Uncertainties in calculated ozone trends need to be re- 
duced. In the 2010 WMO ozone assessment, inconsisten- 
cies in the available ozone profile data prevented a clear 
overall picture of long-term ozone profile changes. A suc- 
cessful assessment of trends for the upcoming 2014 WMO 
ozone assessment requires addressing the problems of time 
series artifacts due to changes in instrumentation, homoge- 
nizing applied satellite retrievals and reducing and describ- 
ing measurement uncertainties. The SPARC/IO3C/IGACO- 
O3/NDACC (SI 2 N) Initiative has brought together scientists 
dealing with ground and satellite ozone profile measurements 
in order to consolidate efforts documenting long-term ozone 
profile changes. SI 2 N has also defined the time frame to en- 
sure that any data set improvements are finalized to ensure 
inclusion in the 2014 WMO ozone assessment analysis. 


This paper provides a summary of the most important 
available ozone data sets, describes measurement methods, 
detection limits and uncertainty estimates, gives details about 
the vertical and temporal resolution of individual data sets 
and provides references to relevant published work. In addi- 
tion, improvements and homogenizations that have occurred 
to date are discussed. This is the first of three overview papers 
discussing the results from the SI 2 N Initiative. The second 
paper describes the validation of the different available ozone 
profile measurements (Lambert et al., 2014) and the third 
paper provides information about and an overview on the 
merging and homogenizing of different data sets and calcula- 
tion of long-term trends (Harris et al., 2014). Ground-based 
ozone profiles, ordered according to the number of available 
stations, are described in Sect. 2. In Sect. 3, satellite-based 
ozone profiles are covered. The satellite measurements are 
organized by measurement techniques; these include solar- 
/stellar-occultation measurements, limb measurements and 
nadir measurements. A summary, including discussion of 
possible temporal and spatial coverage problems for the mea- 
surement systems, is presented in Sect. 4. 

2 Ground-based measurement systems 

The global map with location of the ground-based instru- 
ment systems (ozonesonde, Dobson/Brewer, lidar, FTIR and 
microwave) is shown in Fig. 1. Ozonesondes and Dob- 
son/Brewer instruments provide the longest record of strato- 
spheric ozone variability (see Fig. 2). Ozonesondes are all 
weather sounders, and can measure at night and through 
clouds. Their data quality is good in the lower stratosphere 
and troposphere where satellites measurement quality is less 
reliable, and their vertical resolution is much higher than is 
achievable from remotely sensed measurement. Ozoneson- 
des are inexpensive and easily deployable. Dobson Umkehr 
measurements have captured evolution of the upper strato- 
spheric ozone prior to the satellite era and still remain the 
most inexpensive measurement system requiring minimal 
manpower. Umkehr measurements with Brewer instruments 
were added later, but they offer fully automated and remote- 
controlled operations, and furthermore add spectral UV mea- 
surements. Ozone lidar measurements are self-calibrated and 
provide high vertically resolved measurements with small 
uncertainties (especially for the low to mid-stratosphere). 
Ground-based high-resolution FTIR solar absorption mea- 
surements do not require an absolute calibration either; the 
ozone absorption lines are very narrow and are therefore 
self-calibrated with the reference being the surrounding con- 
tinuum. This technique has also the advantage of provid- 
ing precise total columns in addition to vertical information. 
Ground-based microwave measurements can be made with 
a time resolution of an hour or less throughout a full diur- 
nal cycle; their calibration is based on temperature, the mea- 
surements are insensitive to stratospheric aerosols, and the 
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Figure 1 . World map with locations of long-term ground-based 
measurement sites: ozone sounding sites (red plus signs), lidar sites 
(green stars), microwave sites (light blue triangles), Dobson/Brewer 
sites with Umkehr measurements (yellow squares), and FTIR sites 
(dark blue circles). 

instruments are nearly fully automated and require less per- 
sonnel time and other resources to operate than most other 
ground-based instruments. 

2.1 Ozonesondes 

Ozonesondes are small, lightweight and compact balloon- 
borne instruments, developed for measuring the vertical dis- 
tribution of atmospheric ozone up to an altitude of about BO- 
SS km (see Table la). During flight operation, ozonesondes 
are coupled with standard meteorological radiosondes pro- 
viding temperature, pressure and wind measurements. The 
effective vertical resolution of the ozone profile, which is de- 
pendent on the balloon ascent rate, is about 100-150 m (see 
Table lb and Fig. 3a). 

There were 66 stations contributing data in the 2000s 
(Fig. 1). Ozonesonde records provide the longest time se- 
ries of the vertical ozone distribution throughout both tro- 
posphere and stratosphere; some station records start in 
the 1960s (see Fig. 2 and Table lc). Most of the profiles 
are from the electrochemical concentration cell (ECC)-type 
ozonesonde (Komhyr, 1969) introduced in the early 1970s 
and adopted by most stations in the global network by the 
early 1980s. Almost all data in the most recent decade are 
from ECC sondes. Remaining data are from Brewer-Mast 
(B-M) sondes (Brewer and Milford, 1960), Japanese carbon 
iodine cell sondes (KC96) (Kobayashi and Toyama, 1996), 
and Indian sondes. Prior to the early 1990s, three stations 
in Eastern Europe flew the GDR (German Democratic Re- 
public) sonde. A majority of the data before 1980 is from 
B-M sondes or similar (both the GDR and Indian sondes are 
similar in design to the B-M sonde). All five types use the 
reaction of ozone with potassium iodide in an aqueous solu- 
tion as the method of ozone detection, but the instrumental 


layouts are significantly different (e.g., Smit, 2002; Tarasick 
and Slater, 2008). 

When properly prepared and handled, ECC sondes yield 
profiles with random errors of 3-5 % ( la) and overall uncer- 
tainties of about 5% in the stratosphere (Kerr et al., 1994; 
Smit et al., 2007; Deshler et al., 2008; Liu et al., 2009). Two 
types of ECC sondes are in current use, with minor differ- 
ences in construction and some variation in recommended 
concentrations of the sensing solution and of its phosphate 
buffer. The maximum change in stratospheric response re- 
sulting from these systematic differences is on the order of 
2-3 % (Smit et al., 2007). Other sonde types have somewhat 
larger random errors of 5-10% (Kerr et al., 1994; Parrish 
et al., 2013). 

In early intercomparisons (Attmannspacher and Diitsch, 
1970, 1980), the Indian and the GDR sondes showed sig- 
nificantly larger random errors than other sonde types. The 
largest systematic differences between sondes are in the 
lower stratosphere, where in early intercomparisons the B- 
M and GDR sondes read ~ 5-10 % lower than the ECC and 
KC sondes (Attmannspacher and Diitsch, 1970, 1980). The 
Indian sonde has generally shown little bias in the lower 
stratosphere, but somewhat larger random errors than the 
other non-ECC sondes. In later intercomparisons differences 
were smaller: the KC sondes were biased low by ss 5 % in 
the lower stratosphere, and the B-M sonde generally showed 
a small low bias as well (Kerr et al., 1994; Smit and Kley, 
1998; Deshler et al., 2008). 

In the middle stratosphere (below 28 km) differences in 
sonde response are small. As noted above, both systematic 
and random errors have improved with time, and from 1980 
onwards (between the tropopause and 28 km) the random er- 
ror component of sonde measurements, for all types, is gen- 
erally within ±5 %, and systematic biases between them or 
compared to other ozone sensing techniques are smaller than 
±5 % (SPARC, 1998). 

Above 28 km the measurement behavior of the differ- 
ent sonde types is not consistent due to instrumental un- 
certainties (e.g. pump-flow corrections and sensing solution 
strength) and cannot be generally characterized. Here, B- 
M sondes show increasing underestimation of ozone with 
altitude (De Backer et al., 1998; Stiibi et al., 2008), while 
KC sondes tend to overestimate ozone with increasing alti- 
tude (SPARC, 1998; Smit and Straeter, 2004; Deshler et al., 
2008). In contrast, the performance of ECC sondes between 
28 and 35 km is still good, overestimating ozone compared 
to other measurement techniques by less than 5 % (SPARC, 
1998; Smit et al., 2007). 

However, for ECC sondes recent studies (Johnson et al., 
2002; Smit et al., 2007; Deshler et al., 2008) have demon- 
strated that changes of manufacturer or strength of sens- 
ing solution can introduce significant inhomogeneities (up 
to ±5-10%) in long-term sounding records. Such artifacts, 
like those introduced by changes of sonde type, can be re- 
moved through the use of transfer functions derived from 
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Figure 2. Temporal coverage of the described ground- and satellite-based measurement systems. 


intercomparisons in the laboratory or field. Most station 
long-term records are currently being re-evaluated and ho- 
mogenized under the Ozonesonde Data Quality Assessment 
activity of the SI 2 N Initiative (H. G. J. Smit, personal com- 
munication, 2013). This should both reduce overall uncer- 
tainties and allow increased confidence in trends derived 
from ozonesonde data in future assessments. 

Ozonesonde profiles are archived by the World Ozone 
and Ultraviolet Radiation Data Center (WOUDC), the 
Network for the Detection of Atmospheric Composition 
Change (NDACC) and the Southern Hemisphere Additional 
OZonesondes network fSHADOZ) (see Table Id for more 
details). 

2.2 Dobson/Brewer 

Dobson and Brewer are ground-based spectrophotometer 
instruments. The Umkehr measurement is a sequence of 
morning or afternoon zenith-sky measurements recorded as 
a relative change of transmission at two spectral channels 
(Dobson) or photon counts at individual spectral channels 
(Brewer), all selected in the ultraviolet part of the zenith-sky 
spectrum. Measurements are taken when the sun elevation 
changes between 60° and 90° solar zenith angle (SZA), and 
under the assumption of static atmospheric conditions and 
no clouds in the zenith viewing area. The vertical resolution 
of an Umkehr ozone profile is & 10 km. Ozone profiles (in 
Dobson units (DU)) are historically provided in 10 pressure 


layers, where pressure at the top of the layer is half of the 
pressure at the bottom, while the tropospheric measurement 
is one thick layer below 250 hPa (see Table la, b). The uncer- 
tainty in the retrieved ozone profile is a combination of mea- 
surement error (increases with increasing SZA) and smooth- 
ing error (about 10-20% in the troposphere and low strato- 
sphere and » 5 % in middle and upper stratosphere). 

The Umkehr record length varies by station; the longest 
record is from Arosa, Switzerland (starting in 1956; see Ta- 
ble lc). The latest version of the ozone profile algorithm for 
processing of Dobson Umkehr data (UMK04) is described 
in Petropavlovskikh et al. (2005). Although Umkehr profiles 
show systematic biases in comparisons with other measure- 
ments (Nair et al., 2011; Kramarova et al., 2013b), this is 
of less concern for trend analysis when the UMK04 pro- 
file retrievals are used as monthly mean anomalies. However, 
smoothing errors, especially in the lower stratosphere, create 
low vertical resolution in the Umkehr retrieved profile and 
enhance retrieval noise in the lower stratosphere and tropo- 
sphere. The forward model used to process Brewer Umkehr 
data (03BUmkehr v.2.6, http://www.o3soft.eu/o3bumkehr. 
html) was adapted from the UMK04 model taking into ac- 
count the different optical characteristics of the Brewer in- 
strument (Petropavlovskikh et al., 201 1, see Table Id). 

The Umkehr method uses a technique that minimizes sys- 
tematic errors of instrumental drifts by normalizing to one 
of its measurements. Other methods are used to track instru- 
ment stability, such as once a month checks of wavelength 
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Table la. Summary of measurement techniques and lowest and highest covered level for the five ground-based measurement systems and all 
described satellite measurement systems. 


Name of instrument 

Measurement technique 

Covered vertical region: lowest level 

Covered vertical region: highest level 



Ground-based measurement systems 


Ozonesondes 

Electrochemical conversion of ozone 

Surface 

30-35 km altitude (bursting point of balloon) 

Dobson/Brewer 

Umkehr, ground-based, UV part of the spectrum 

Tropospheric layer: surface to ~ 250 hPa 

To top of atmosphere 

Lidar 

Differential Absorption Lidar (DIAL), ground-based, 
UV laser radiation 

Tropopause 

45-55 km depending on the lidar systems 

FTIR 

Ground-based solar absorption measurements in the infrared. 
The ozone profile is retrieved from the pressure-broadened line 
shape. 

Surface 

40-45 km 

Microwave 

Measurement of the emission spectrum of a thermally ex- 
cited rotational transition of ozone in the millimeter wave- 
length range. The ozone profile is retrieved from the pressure- 
broadened line shape. 

Typically Rs 20 km 

Varies by spectral line frequency. Up to 72 km at night, and 
66 km during the day for the transition at 1 10.836 GHz. 


Satellite-based measurement systems 


SAGE 

Solar occultation 

Tropopause or 10 km whichever is lower 

55 km 

HALOE 

Solar occultation instrument; 9.4 to 10.4 pm ozone channel 
(SPARC, 1998) 

Cloud top or to just below the tropopause 

To near 0.005 hPa for temperature and ozone 

ILAS/ILAS-II 

Satellite-based, solar occultation 

Cloud top: 10 km 

To top of atmosphere 

POAM 

Solar occultation in the UV-V is spectral range 

5 km (or cloud top) 

60 km 

ACE-MAESTRO 

Solar occultation spectrophotometry in the Chappuis band 

Rs 5 km 

35 km 

ACE-FTS 

Satellite-based, solar occultation, infrared spectrum 

Cloud tops (~ 5 km) 

95 km 

GOMOS 

Satellite-based stellar occultation instrument in UV-VIS-NIR 

Typically cloud top, however, lowest valid altitude ~ 15 km 

100 km 

LIMS 

Satellite-based instrument, limb-infrared emission experiment; 
9 to 10 micrometer band for ozone 

250 to 300 hPa or higher cloud top 

To near 0.01 hPa for T and ozone 

UARS MLS 

Microwave limb emission sounding (satellite) 

100 hPa 

0.2 hPa 

Aura MLS 

Microwave limb emission sounding (satellite) 

215 hPa 

0.02 hPa 

MIPAS 

Satellite-based mid-infrared limb emission 

It depends on the measurement mode. Most data are in nominal 
observation mode: 

6 km (or cloud top) for 2002 to 2004 

7-12 km (following tropopause height) or cloud top for 

2005 to 2012; 

Middle/upper atmosphere modes: above 20 or above 40 km 
These are nominal numbers; real numbers can vary. 

« 68/72 km (2002-2004/2005-2012) in nominal observation 
mode; up to 100 km in middle/upper atmosphere modes. 

Odin-SMR 

Microwave limb emission sounding (satellite) 

*=» 12 km (poles), 17-1 8 km (mid-latitudes) 

Rs 60 km 

SMILES 

Submillimeter-wave limb emission sounding (onboard 
the International Space Station) 

16km 

Rs 95 km 

SCIAMACHY 

UV-V is limb scatter 

Tropopause ( 1 0 km) 

Mesosphere (60 km) 

Odin-OSIRIS 

Satellite based, UV-visible, limb scatter 

Cloud tops 

55 km 

SBUV 

Satellite-based nadir-view backscattered UV radiation 

Surface 

Top of atmosphere 

GOME/GOME-2 

Optimal Estimation method, satellite-based instrument looking 
in nadir direction, UV-VIS part of the spectrum 

Surface 

Top of atmosphere 

OMI 

UV-VIS nadir spectrometer 

Surface 

Up to 0.3 hPa 

IASI 

Fourier transform spectrometer, thermal infrared 

Tropospheric layer or lower if thermal contract is good 

Total columns (sensitivity up to around 45 km) 


registration by using standard discharge lamps, ratio of mea- 
surements at several spectral channels (for Brewer instru- 
ments, daily), and optical wedge calibration produced dur- 
ing Dobson intercomparisons (every 4 years). These meth- 
ods track degradation of the optical system. Effects of ozone 
cross-section values and uncertainties of stray light effects 
on Dobson and Brewer Umkehr ozone profiles are discussed 
in Petropavlovskikh et al. (2011) and WMO (2008). The re- 
placement of an instrument at the station can cause a step 
change in the Umkehr record due to the difference in the 
out-of-band contribution that is specific to each instrument. 
Thus, additional homogenization of the record may be re- 
quired (Zanis et al., 2006; Miyagawa et al., 2009). In addi- 
tion, interference from optically thick aerosol loading (vol- 
canic) at Umkehr spectral channels results in systematic er- 
rors that can be as large as 10-15 %. Therefore, about two 
years of data are typically removed after the Pinatubo (1991) 
and El Chichon (1982) volcanic eruptions. Data from 28 
stations with current Umkehr measurements and of 58 sta- 
tions with historical Umkehr measurements are archived at 
WOUDC, Toronto, Canada (see Table Id for the URL). 


2.3 Lidar 

Lidar (light detection and ranging) is an active, remote- 
sensing instrument that uses the interaction between a laser 
beam and atmospheric molecules and particles. Ozone lidar 
measurements are performed using the Differential Absorp- 
tion Lidar method (DIAL), which is based on the absorption 
of ultraviolet radiation by ozone and requires the emission 
of two laser wavelengths (so-called “on” and “off" wave- 
lengths) characterized by different ozone absorption cross 
sections (Pelon et al„ 1986; Godin et ah, 1989; McDermid 
et ah, 1990). The use of pulsed lasers provides range re- 
solved measurements. Ozone DIAL systems include one or 
two lasers, depending on the technique used for the gen- 
eration of the off wavelength, an optical receiving system 
that includes a telescope for the collection of the laser light, 
a spectrometer or filters for the separation of the received 
wavelengths and an electronic data acquisition system. In 
the case of stratospheric ozone, photon counting is generally 
used for the acquisition of the lidar signal in order to provide 
high sensitivity and low noise. The ozone number density is 
commonly retrieved from the difference of the slopes of the 
logarithm of the lidar signals corrected for the background 
noise and saturation affecting large signals originating from 
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Table lb. Summary of the vertical resolution, native vertical grid and native ozone units for the five ground-based measurement systems and 
all described satellite measurement systems. 


Name of instrument 

Vertical resolution* 

Representation grid 

Native vertical grid* 

Native ozone units* 



Ground-based measurement systems 


Ozonesondes 

100-150 m 

(at ascent rate of 5 m s - 1 and response time of 
20-30 s) 

1 00-200 m 

Flight time 

Partial ozone pressure in milli-pascal (mPa) 

Dobson/Brewer 

«= 10km 

R3 10 km 

Pressure levels 

Dobson units (DU) 

Lidar 

The resolution varies as a function of altitude 
from 0.5 km below 20 km, 2 km around 30 km 
to more than 5 km above 45 km. 

Depending on the systems: 0.15 or 0.3 km 

Altitude levels 

Ozone number density (cm -3 ) 

FTIR 

8-10 km for the ground ~ 9 km layer 

353 8 km for the lower and middle stratospheric 
layer 

s» 15-20 km for the upper stratospheric layer 
(S3 28-45 km) 

In NDACC HDF archived files: the retrieval 
grid contains, depending on the station, about 
41 to 47 levels (from ground to 100 km). 

Pressure levels 

Volume mixing ratio (ppmv) 

(provided also as partial columns (mol cm -2 )) 

Microwave 

Varies by instrument, typ. 8-1 0 km, from 
«3 20 to 40 km, then increasing to ss 17 km at 
«3 60 km, see Sec. 2.5 

Varies by instrument. Typ. «3 2 km intervals 

Pressure levels 

(grid on which data are provided varies by 

Volume mixing ratio (ppmv) 




instrument) 




Satellite-based measurement systems 


SAGE 

1km 

1 km 

Altitude levels 

Ozone number density (cm -3 ) 

HALOE 

Measurement: 1.8 km instantaneous field of 
view (IFOV); retrieval resolution is 2.3 km 

2.3 km vertical resolution 

Solar zenith angle or apparent altitude, sampled 
every 0.3 km. 

(provided also on pressure and associated 
hydrostatic altitude levels) 

Transmission vs. solar zenith angle or apparent 
altitude 

(provided as volume mixing ratio (ppmv)) 

ILAS/ILAS-II 

ILAS: % 1.9 km 
ILAS-II: « 1.6 km 

Rsl km 

Tangent height (in km) 

Volume mixing ratio (ppmv) 

POAM 

ss 1.5 km 

«3 1 km 

Altitude levels 

Ozone number density (cm -3 ) 

ACE-MAESTRO 

S3 1.2 km 

R3 1.0 km 

Tangent levels 

Concentration (mol cm -2 ) 

(provided as volume mixing ratio (ppmv)) 

ACE-FTS 

ss 3—4 km 

1 km (interpolated to standard grid) and R3 3- 
4 km (as measured) 

Altitude levels 

Volume mixing ratio (ppmv) 

GOMOS 

Resolution of ozone measurements: 

2 km below 30 km 

3 km below 40 km 

smooth transition between 30-40 km 

Vertical sampling resolution: about 0.5-1. 7 km. 
Variable for each profile. 

Altitude levels 

Ozone number density (cm -3 ) 

LIMS 

Measurement: 1.8 km instantaneous field of 
view (IFOV); Retrieval: 3.7 km 

Level 2 data are provided every 0.88 km, having 
a vertical resolution of ^ 3.7 km 

Pressure levels 

(data are also provided on altitude levels) 

Volume mixing ratio (ppmv) 

UARS MLS 

3.5— 4km in the stratosphere and 5-8 km in the 
mesosphere 

3 km (6 levels per decade change in pressure for 
stratosphere) 

Pressure levels 

Volume mixing ratio (ppmv) 

Aura MLS 

2.5-3 km in upper troposphere to lower meso- 
sphere, 4 to 6 km in upper mesosphere 

R3 3 km (6 levels per decade change in pressure) 
for pressures larger than 0. 1 hPa; grid resolution 
degrades by a factor of two for pressures less 
than 0.1 hPa. 

Pressure levels 

Volume mixing ratio (ppmv) 

MIPAS 

Optimized resolution nominal mode (years 
2005-2012) 

IMK: 2.4— 3.6 km (reduced resolution nominal 
mode) 

ESA processor: 2 km at 10 km, 3.2 km at 30 km, 
5 km at 70 km, mostly independent on atmo- 
spheric conditions 

Full resolution mode (years 2002-2004) 

IMK: 3.5km at 10km, 5km at 42km, 8km at 
60 km. 

ESA: 4 km at 10 km, 5 km at 30 km, 8 km at 
70 km. 

IMK: retrieval and data representation grid for 
reduced resolution nominal mode is 1 km from 
4—44 km, 2 km from 44-70 km; 5 km above 
ESA processor: grid of the measured tangent 
altitudes 

MIPAS2D: 3 km from 6 to 42 then 47, 52, 60 
and 68 km 

IMK: Altitude levels 

(ESA: provides pressure levels where the re- 
trieved pressure corresponds to the measured 
tangent altitude) 

Volume mixing ratio (ppmv) 

Odin-SMR 

2.5-3.5 km 

1 .5 km and 3 km 

Pointing information is retrieved from pressure 
broadened spectral lines 
(provided on altitude levels: altitude-pressure- 
temperature information provided for interpola- 
tion between grids) 

Volume mixing ratio (ppmv) 

SMILES 

2-4 km 

2.5 km 

Altitude levels 

Volume mixing ratio (ppmv) 

SCIAMACHY 

3-5 km 

1 km 

Geometric height (km) 

Ozone number density (cm -3 ) 

Odin-OSIRIS 

*3 1 .5 km in UTLS and 2 km higher up 

3« 1 km 

Altitude levels 

Ozone number density (cm -3 ) 

SBUV 

*3 6 km at 3 hPa 

* 3 15 km in troposphere and lower mesosphere 

3.2 km 

Pressure levels 

Dobson units (DU) 

(provided also as volume mixing ratio (ppmv)) 

GOME/GOME-2 

*3 7-15 km 

353 4 km 

Pressure levels 

Dobson units (DU) 

OMI 

6-10 km 

2-5 km 

Pressure levels 

Dobson units (DU) 

IASI 

«3 7 km 

1 km up to 40 km 

2 km from 40 km to 60 km 

Altitude levels 

Ozone number density (cm -3 ) 

(provided as volume mixing ratio (ppmv) or 





ozone partial columns (DU)) 


* specified by instrument/averaging kernel 
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Table lc. Summary of the covered time period, the number of measurements and the time for measurements for the five ground-based 
measurement systems and all described satellite measurement systems. 


Name of instrument 

Covered time period 

Average number of measurements 
(per day/per week/per month) 

Measurement time of day 

Ground-based measurement systems 

Ozonesondes 

Varies from station to station, longest record from mid 1960s, 

1-2 vertical ozone soundings per week 

Ozonesondes are usually launched at daytime before about 


early 1970s 


noon time 

Dobson/Brewer 

Varies from station to station, longest record from Arosa, 

Average is about 15 per months (~ 200 per year; depending on 

Sunrise and sunset between 70 and 90° solar zenith angle (SZA) 


Switzerland: 1956-present 

station and season) 

(only with clear sky) 

Lidar 

Varies from station to station, the longest records start in the 

Up to 20 measurements per month (depending on weather con- 

Nighttime (only with clear sky) 


late 1980s (OHP, France; Hohenpeissenberg, Germany; Table 
Mountain, USA) 

ditions) 


FTIR 

Varies from station to station, usually 1995-present. 

Mean of all stations: ~ 2.5 per day (7 per week; 15 per month) 

Daytime (only with clear sky). 

Microwave 

Sep 1992-present (Lauder), Jun 1995-present (Mauna Loa), 
2000-present (Payeme). Other records back to 1990’s may be 
available. 

4 to 48 profiles per day in suitable weather, depending on station 

Day- and nighttime (only with clear sky) 

Satellite-based measurement systems 

SAGE 

SAGE I: Feb 1979-Nov 1981 
SAGE II: Oct 1984-Jul 2005 
SAGE III: Dec 2001-Mar 2006 

Up to 15 sunrise and 15 sunset events per day 

Sunrise/sunset 

HALOE 

Oct 1991-Nov 2005 

Up to 15 sunrise and 15 sunset events per day 

Sunrise/sunset 

ILAS/ILAS-II 

ILAS: 30 Oct 1996-29 Jun 1997 

ILAS: 28 profiles per day (14 sunrise and 14 sunset profiles) 

Sunrise/sunset 


ILAS-II: 2 Apr 2003-24 Oct 2003 

ILAS-II: 28 profiles per day (14 sunrise and 14 sunset profiles) 


POAM 

POAM II: Sep 1993-Nov 1996 

14 measurements per day in each hemisphere 

Measurement performed at local sunset in NH, but in the SH 
they are at local sunset in the winter and local sunrise in the 




summer. 


POAM III: Apr 1998-Nov 2005 



ACE-MAESTRO 

Feb 2004-present 

Up to 30 profiles per day 

Sunrise/sunset 

ACE-FTS 

Feb 2004-present 

Up to 30 profiles per day 

Sunrise/sunset 

GOMOS 

Aug 2002-Apr 2012 

Daily: 100-200 profiles 

Nighttime (around 22:00) 


Reduced number of measurements due to 

Monthly: 3000-6000 profiles 



instrumental problems in year 2005 and since Oct 2011. 

Yearly: 36000-72000 profiles 


LIMS 

Oct 1978-May 1979. 

rs 3000 profiles per day (the LIMS observation duty cycle was 

LIMS operated in a sun-synchronous mode, taking its tangent 



of the order of 5 days on and one day off for most of its mission 

point scans at about 1 300 and 2300 h (at the Equator), but with 



life) 

little change of local time with latitude until reaching high lati- 
tudes. 

UARS MLS 

Sep 1991-Jul 1999 

R» 1350 profiles/day 

Day and night measurements are routinely performed around 
the orbit, with all local times covered at any given latitude dur- 
ing each satellite yaw cycle (every 36 days on average) 

Aura MLS 

Aug 2004-present 

rs 3500 profiles/day 

Two local times are sampled every orbit at any given latitude, 
daytime and nighttime, with local times typically near 13:00- 
14:00 and near 01:00-02:00 for non-polar latitudes 

MIPAS 

Jun 2002-Apr 2012 with gaps in 2004 

Rs 1400 profiles/day (optimized resolution phase, 
nominal mode) 

rs 1050 profiles/day (full resolution phase, nominal mode) 

10:00 and 22:00 local time 

Odin-SMR 

Nov 2001 -present 

In the stratospheric mode, Odin measures up to (almost) 1000 
vertical profiles per observation day based on up to rs 65 scans 
per orbit and 14—15 orbits per observation day. There were 10 
observation days per month before Apr 2007 (every third day), 
15 days per month since 2007 (every other day) 

Day- and nighttime (thermal emission measurement). 

SMILES 

Oct 2009-Apr 2010 

Rs 1600 per day 

Day- and nighttime measurements 

SCIAMACHY 

Aug 2002-Apr 2012 

rs 1400 measurements per day 

Daytime measurements only. 

Odin-OSIRIS 

Nov 2001 -present 

OSIRIS duty cycle on Odin changed in 2007, so scans (between 

Only a.m. data are included in the SI 2 N analysis, typically mea- 



—60 and 60°) are reported in pairs (pre-2007 and post-2007): 
Average of 75/150 per day (450/900 per week; 2000/4000 per 
month) 

sured between 05:30 and 07:30, local time. 

SBUV 

Oct 1978-present continuous coverage; 1970-1976 NOAA-4 
coverage, partial coverage after mid-1972. 

Rs 1100 per day 

Daytime measurements only. 

GOME/GOME-2 

GOME: 1995-2003 
GOME-2A: 2006-present 
GOME-2B: 2012-present 

Rs 20 000 measurements per day 

Daytime measurements only. 

OMI 

2004-present 

rs 140000 measurements per day 

Daytime measurements only. 

IASI 

2007-present 

rs 1.2 x 10 6 measurement per day (about 50 % are cloud 

Two overpasses per day at every location (around 09:30 and 



contaminated) 

21:00, local time) 


lowermost altitudes. A correction is applied corresponding to 
differential terms linked to molecular and aerosol scattering, 
and absorption by other atmospheric constituents (e.g. NO 2 , 
S0 2 ). 

The DIAL technique does not require any calibration and 
the laser wavelengths are chosen such that the correction 
term represents less than 10% of the ozone number den- 
sity derived from the slope of the lidar signal in the alti- 
tude range of interest. When stratospheric aerosol loading 
is large, such as after cataclysmic volcanic eruptions, this 
condition is not fulfilled and DIAL measurements are locally 
perturbed at the altitude of the aerosols. In such conditions, 
two additional wavelengths corresponding to the first Stokes 


vibrational Raman scattering of the laser radiation by atmo- 
spheric nitrogen can be detected and the Raman lidar signal 
is used to retrieve ozone (McGee et al., 1993). 

Most lidar measurements are performed during the night 
and averaged over several hours, resulting in a horizontal 
resolution of 50 to 250 km over the whole altitude range 
of the measurement, depending on atmospheric conditions. 
The vertical resolution decreases as a function of altitude 
ranging from several hundred meters at lower altitudes to 
several kilometers above 40km (Table la, b; Fig. 3a). The 
systematic uncertainty ranges from a few percent below 
20 km to more than 10-15 % above 45 km (Godin-Beekmann 
et al., 2003). Ozone lidar measurements are self-calibrated. 
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Table Id. Summary of the latest data version, URLs where data are available, and information about references in the S1~N special issue for 
the five ground-based measurement systems and all described satellite measurement systems. 


Name of instrument 

Latest version of data 


Where data are available 

Reference in SI 2 N special issue 

Ground-based measurement systems 

Ozonesondes 

Outcome of 03S-DQA activity as part of SI 2 N 


WOUDC, Toronto, Canada: 

http://www.woudc.org/data_e.html 

NDACC website (selected stations): 

http://www.ndacc.org 

SHADOZ website (tropical stations): 

http://croc.gsfc.nasa.gov/shadoz/ 

(not available at this point) 

Dobson/Brewer 

Described in Petropavlovskikh et al. (2005) 


WOUDC, Toronto, Canada: 
http://www.woudc.org/data_e.html 

“ 

Lidar 

Depends on each lidar system. 


NDACC website: 
http://www.ndacc.org 

Nairet al. (2013) 

FTIR 

Described in Vigouroux et al. (2008) 


NDACC website: 
http://www.ndacc.org 

(not available at this point) 

Microwave 

Version 6. 


NDACC website: 
http://www.ndacc.org 

Special Lauder and Mauna Loa hourly processed data sets available on request 
from PI. 

Smit and Kley (1998) 

Satellite-based measurement systems 

SAGE 

Version 7 


Atmospheric Science Data Center, NASA Langley. 
http://eosweb.larc.nasa.gov/project/sage2/sage2_v7_table 

Damadeo et al. (2013) 

HALOE 

Version 19 ozone is described briefly in Randall et al. (2003). 


NASA Goddard Earth Sciences Data & Information Services Center (GES 
DISC) : http://daac.gsfc . nasa. gov/ 



More details: Briihl et al. (1996) and SPAPC (1998). 


HALOE Webpage: http://haloe.gats-inc.com 


ILAS/ILAS-II 

ILAS: Version 6.10 


ILAS : http://db.cger.nies.go.jp/ilas/ 

- 


ILAS-II: Version 3.0 


ILAS-II: http://db.cger.nies.go.jp/ilas2/ 


POAM 

Latest retrieval version is 4. Version 4 ozone differs very little from version 

NASA Langley Atmospheric Sciences Data Center: http://eosweb.larc.nasa. 

- 


3. The version 3 algorithm and error analysis are described in 

Lumpe et al. 

gov/project/poam3/poam3_table 



(2002). Changes for the version 4 algorithm are described in a pdf document on 




the POAM ftp archive. 


Data also available through Naval Research Laboratory: 
and: http://www.nrl.navy.mil/rsd/7220/poam 


ACE-MAESTRO 

Described in McElroy et al. (2007) 


ACE-MAESTRO: From ACE team, University of Waterloo, Waterloo, Canada 
(for most recent version) 

Publicly available data sets from: 
http://www.ace.uwaterloo.ca/data.html 


ACE-FTS 

Described in Boone et al. (2005) 


ACE-FTS: From ACE team. University of Waterloo, Waterloo, Canada (for v3.5 
- most recent versions) 

(v2.2 + updates - Publically available data sets from http://www.ace.uwaterloo. 
ca/data.html) 

(not available yet) 

GOMOS 

GOMOS IPF 6.01 released in 2012 


ESA EARTHNET: 

(Kyrola et al., 2013) 




https://earth.esa.int/web/guest/missions/esa-operational-eo-missions/envisat/ 

instruments/gomos 

Vertically gridded data set: 





(Sofieva et al., 2013) 

LIMS 

Described in Remsberg et al. (2007) 


NASA Goddard Earth Sciences Data & Information Services Center (GES 
DISC): 

http://daac . gsfc .nasa.gov/ 


UARS MLS 

Version 5 was the last major public release of UARS MLS O3 data (Livesey 

The official public distribution location for UARS MLS data is at the NASA 

(not available yet) 


et al., 2003) 


GES-DISC Mirador website: 
http://mirador.gsfc . nasa. gov 


Aura MLS 

v2.2 


The official public distribution location for Aura MLS data is at the NASA 

GES-DISC Mirador website: 

http://mirador.gsfc.nasa.gov 

(not available yet) 

MIPAS 

IMK: V5_03_221 (von Clarmann et al., 2009; Laeng et al., 2014) 

IMK, Germany: 

Eckert et al. (2014) 


ESA: ML2PP V6 


http://www.imk-asf.kit.edu/english/308.php 



MIPAS2D - Version 2.3 (see Dinelli et al., 2010, 2012) 


Oxford University, UK: 
http://www.atm.ox.ac.uk/group/mipas/L20XF/ 

ESA: 

https://earth.esa.int/web/guest/data-access/browse-data-products 
MIPAS2D - ISAC-CNR, Italy: 
www.isac.cnr.it/~rss/mipas2d.htm 

Raspollini et al. (2013b) 

Odin-SMR 

Version 2.1. 


http://odin.rss.chalmers.se 

- 

SMILES 

JAXA: 


JAXA: 

JAXA: 


Version 3.0 (118-12-0602) 


ISAS/JAXA DARTS 

(not available yet) 


Described in Mitsuda et al. (201 1) and Imai et al. (2013) 


http://darts.isas.jaxa.jp/iss/smiles/ 



NICT: 


NICT: 

NICT: 


Version 2. 1 .5 


https://smiles-p6.nict.go.jp/products/research_latitude-longitude.jsf 

Kasai et al. (2013) 


Described in Baron et al. (201 1); Kasai et al. (2013) and Kreyling et al. (2013) 



SCIAMACHY 

Version 2.5 

Described in Rahpoe et al. (2013) and Mieruch et al. (2012) 


http://www.iup.uni-bremen.de/ 

Rahpoe et al. (2013) 

Odin-OSIRIS 

Described in Degenstein et al. (2009) 


http://odin-osiris.usask.ca/?q=node/280 

(not available yet) 

SBUV 

Version 8.6 


Goddard DISC 

Bhartia et al. (2013); DeLand 
et al. (2012); Kramarova et al. 
(2013a, b) 


Described in Bhartia et al. (2013) and in McPeters et al. (2013) 


http://disc.sci.gsfc.nasa.gov/daac-bin/DataHoldingsMEASURES.pl? 





PROGRAM_List=RichardMcPeters 


GOME/GOME-2 



GOME-2: EUMETSAT archive: 
http://archive.eumetsat.int/ 

GOME-2: TEMIS: 
http://www.temis.nl/profiles/ 

All data: In the future via Ozone_cci web portal (ESA): 
http://www.esa- ozone- cci .org 


OMI 

Version 1.1.1 


OMI: Goddard DISC 
http://disc.sci.gsfc.nasa.gov/ 

- 

IASI 

Described in Hurtmans et al. (2012) 


Upon request. Will be operationally distributed by Eumetsat in 2014 

- 


Described in Du four et al. (2012) 


Eumetsat data archive: 





www.eumetsat.int 
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Figure 3. Vertical resolution shown as function of altitude for each of the described instruments systems: (a) ground-based instruments, 
(b) solar/stellar occultation instruments, (c) nadir-measuring instruments, and (d) limb-measuring instruments. The SMILES profile in panel 
(d) is shown as a brown dashed line that is plotted on top of the ODIN-SMR (light purple). Note that there are two profiles shown for 
MIPAS in panel (d): the solid purple line represents the KIT retrieval, and the dash-dotted purple line represents the Bologna retrieval. The 
oscillating altitude resolution of the MIPAS KIT profile is a consequence of the evaluation of the averaging kernels on a vertical grid finer 
than the tangent altitude spacing: the vertical resolution is best at the tangent altitudes and coarsest between the tangent altitudes. Retrieval 
resolution for the Oxford and ESA processors is comparable to the KIT results; however, a courser grid is used for the retrieval and thus it 
does not show the oscillating features of the KIT processor. For more detail on the different MIPAS retrievals see Sect. 3.3.4. 


However, measurement stability can be affected by instru- 
mental artifacts such as alignment error, changes in the inten- 
sity of the lidar signals due to varying laser power, or changes 
in weather conditions (Godin-Beekmann et al., 2003). 

There are twelve stratospheric ozone lidar instruments op- 
erational globally at this point (see Fig. 1), with the longest 
continuous records starting in the late 1980s (Steinbrecht 
et al., 2009). Data are archived in the NDACC database (see 
Table Id for the URL). 

2.4 FTIR 

Ground-based FTIR (Fourier transform infrared spec- 
troscopy) solar absorption measurements are performed over 
the 600^-500 cm -1 spectral range, under clear-sky condi- 
tions using primarily high-resolution spectrometers such as 
the Broker 120M (or 125M) or Broker 120HR (or 125HR). 
The spectrometers can achieve a spectral resolution of 0.0035 
and 0.002 cm -1 , respectively. The advantage of the FTIR 
technique is that, for atmospheric gases such as ozone which 
have very narrow lines, an absolute calibration is not needed; 
the ozone absorption signatures are self-calibrated with the 
reference being the surrounding continuum. Nevertheless, 


for the determination of absorber profiles and of vertical 
trends, especially in the upper-stratosphere, it is important 
to regularly monitor the instrumental line-shape. This is 
achieved with gas cell measurements (Hase et al., 1999). 

In addition to total columns, low-vertical-resolution pro- 
files can be obtained from the temperature and pressure de- 
pendence of the absorption line shapes. Two profile retrieval 
algorithms are widely used, PROFITT9 (Hase, 2000) and 
SFIT2 (Pougatchev et al., 1995); both are based on the op- 
timal estimation method (Rodgers, 2000). A European FTIR 
network has been optimized by using a common retrieval 
strategy to derive ozone trends (1995-2005) in 4 layers 
(Vigouroux et al., 2008), one in the troposphere and three 
in the stratosphere, up to about 45km (Table la, b). Thus, 
a high degree of freedom for the signal (DOFS ~ 4.5) is 
achieved with the use of the 1000-1005 cm -1 window (Bar- 
ret et al., 2002). The random error on total columns is about 
2^1 % (Garcia et al., 2012; Vigouroux et al., 2008); the lead- 
ing error source being the temperature. The estimated ran- 
dom error on the ozone partial columns can vary a little de- 
pending on the station (different instrument, ozone natural 
variability impacting the smoothing error). For example, at 
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Kiruna the random error for the layers ground- 10 km, 10- 
18 km, 18-27 km, and 27^-2 km are 10, 7, 9, and 6%, re- 
spectively (Vigouroux et al., 2008). The smoothing error is 
dominant in the two lowest layers, whereas the temperature 
error dominates in the middle and upper stratosphere. Further 
improvement (on-going research, not applied to the NDACC 
network yet) can be achieved by retrieving a simultaneous 
temperature profile, thus reducing the total column random 
error to less than 1 %, and the upper stratospheric error to 
3 % (Garcia et al., 2012). The systematic error is dominated 
by the uncertainties of the spectroscopic parameters (line in- 
tensity and air broadening coefficient). The systematic error 
due to the line intensity uncertainty was estimated to 2-3 % 
(Barret et al., 2003; Schneider et al., 2008) on total and par- 
tial columns. The error due to the air broadening coefficient 
uncertainty is negligible for total columns but can be as large 
as 4% when partial columns are concerned (Barret et al., 
2003; Schneider et al., 2008). Other sources of systematic 
errors, such as the temperature or the instrumental line-shape 
(if not correctly determined), should also be considered in 
the upper layers (Garcia et al., 2012). 

The ozone trends in Vigouroux et al. (2008) were calcu- 
lated with a bootstrap re-sampling method, applied to the 
daily mean partial columns of the ozone time series. This 
work was extended up to 2009 in the WMO 2010 ozone 
assessment report (WMO, 2011) for five European stations 
(28° N to 79° N). A current update up to 2012 is ongoing, in- 
volving additional stations, including stations in the Southern 
Hemisphere. The spectroscopic database has been changed 
to HITRAN 2008 (Rothman et al., 2009), and an additional 
effort of homogenization has been made in contrast to the 
work described in Vigouroux et al. (2008); all stations are 
using a priori information from the model WACCM. 

The length of the FTIR ozone time series varies by sta- 
tion (see Fig. 1 for the geographical distribution of the sta- 
tions); the longest time series available starts in 1995, and 
the shortest in 2002 (Fig. 2; Table lc). In the coming years, 
more recently established observatory stations (e. g. Reunion 
Island: measurements started in 2009; Toronto: regular data 
containing the 1000-1005 cm -1 range started in 2010) or sta- 
tions that have not yet reprocessed their retrievals using the 
optimized common strategy (e.g.: Eureka, Bremen) could be 
included in ozone trends studies. 

The FTIR data not yet available on the NDACC database 
(see Table Id for the URF) can be acquired by direct con- 
tact with the responsible principal investigator. There are cur- 
rently 17 active long-term NDACC FTIR sites. An additional 
long-term site (Kitt Peak, AZ, USA (31.9° N, 111.6°W)) 
has data from 1978-2005, but operations have unfortunately 
ceased. 

2.5 Microwave 

The NDACC microwave ozone profiling instruments mea- 
sure the spectra of emission lines produced by thermally 


excited, purely rotational ozone transitions at millimeter 
wavelengths. All are based on the same principles but dif- 
fer in technical details. All use a sensitive heterodyne down- 
converter (receiver) that produces a replica of the spectrum 
it receives from the sky at a much lower so-called interme- 
diate frequency where it can be processed by a filter bank 
or digital FFT spectrometer. In all systems, the spectral in- 
tensity scale is established by substituting the thermal radia- 
tion from two black body sources for the radiation from the 
sky at the receiver input. One source is at ambient tempera- 
ture, the second is typically chilled using liquid nitrogen. The 
ozone altitude distribution is retrieved from the details of the 
pressure broadened line shape, typically using various im- 
plementations of the optimal estimation method of Rodgers 
(1976). The attenuation of the ozone signal in the tropo- 
sphere is determined by measuring the tropospheric thermal 
emission and relating the tropospheric opacity to its emis- 
sion using a radiative transfer model. Stratospheric aerosols 
do not affect microwave measurements because the wave- 
length is large compared to the aerosol particle size. The 
altitude range is between about 20 and 55-72 km, depend- 
ing on the instrument (Table la, b). The fundamental na- 
tive measurement units are mixing ratio vs. pressure. The 
instruments operate continuously, and profile retrievals are 
obtained in weather ranging from clear to some overcast con- 
ditions. Temporal resolution can be hourly or less, so diurnal 
variations of ozone in the stratosphere and mesosphere can 
be observed. 

The instruments at Fauder (New Zealand), Bern and Pay- 
erne (Switzerland), and Mauna Foa/Hawaii (USA) have 
operated essentially continuously from 1992, 1994, 2000, 
and 1995 to the present respectively. Ny Alesund (Nor- 
way) has been making measurements since 1994, with some 
extended breaks. These are NDACC Primary Instruments 
listed in http://www.iapmw.unibe.ch/research/collaboration/ 
ndsc-microwave/instruments/index_spe.html. Other instru- 
ments, also listed on the above site, have operated 
over periods of months to a decade (Table lc). Instru- 
ment and measurement details for Fauder and Mauna 
Foa were published in Parrish et al. (1992), for Pay- 
erne in Maillard-Barras et al. (2009) and Hocke et al. 
(2007), for Bern in Studeret al. (2013), for Ny Alesund in 
Palm et al. (2010). Data from these instruments are archived 
on the NDACC webpage (see Table Id for URF). 

Theoretical investigations of microwave ozone measure- 
ment errors and vertical resolution have been reported by 
Connor et al. (1995) and Palm et al. (2010). The former 
is applicable to the Fauder and Mauna Foa measurements, 
the latter to the Ny Alesund measurements. As the ver- 
tical resolution of microwave measurements can vary by 
amounts approaching a factor of two from one instrument 
to another, it is important to download the averaging kernels 
along with the data and consider them in the analysis when 
comparing microwave measurement results. Discussions of 
the averaging kernels and vertical resolution are found in 
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Connor et al. (1995) and Caliesi (2000). The values shown 
in Fig. 3a are an ensemble average. Because instrument cal- 
ibration is based on temperature it can theoretically be kept 
stable over long periods of time. In practice, a more impor- 
tant issue is that the details of the spectral line shape must 
be measured with a precision that is technically challenging 
to achieve. The types of specUal errors most commonly en- 
countered tend to propagate into profile errors in the lower 
to middle stratosphere. Despite these challenges, the pro- 
files of the Mauna Loa and Lauder instruments were found 
to agree with others to within < 5 % from & 22 to 65 km 
during formal blind intercomparison campaigns involving 
several types of instruments (McDermid et al., 1998a, b; 
McPeters et al., 1999). Differences between the Bern (Studer 
et al., 2013) instrument and others mostly fall within 7 %, 
as do those at Payerne up to ~ 50 km. Relative drifts be- 
tween the Lauder and MLO measurements and those made 
with lidar, ozonesonde, SAGE II, and HALOE were typi- 
cally < 0.5 %yr -1 from 22 to 60 km over the period from 
the mid-1990s to the mid-2000s (Boyd et al., 2007). While 
no decadal-scale drift tests have been reported for the other 
stations, Steinbrecht et al. (2009) found that the Bern mea- 
surements were very consistent with others over the period 
from 1994 to 2008. 


3 Satellite-based measurement systems 

In the following paragraphs several satellite-based systems 
measuring ozone profiles are described. Systems are grouped 
according to measurement technique. The solar occultation 
technique provides specUal measurements with a very high 
signal-to-noise ratio allowing the detection of species with 
low concenUations without the need for averaging measure- 
ments. This self-calibrating method is less susceptible to 
changes in instrument performance over the mission life- 
time because of the reference exo-atmospheric measure- 
ments used in the retrieval. Stellar occultation measurements 
are relative by definition and therefore they are less sensitive 
to instrumental degradation. Due to the multitude of stars, 
global coverage is obtained. Since stars are point sources, 
the pointing information is excellent (uncertainty < 30 m) 
and the vertical resolution very good, about 2 km around the 
ozone maximum and below. The limb viewing technique of- 
fers high sensitivity to a number of trace gases with vertical 
resolution between 1 and 3 km (depending on instrument and 
wavelength) and provides day and night retrievals (only day 
retrievals for limb-scatter instruments), with dense sampling 
along the sub-orbital ttacks and hundreds to a few thousands 
of profiles per day. Self-calibration is typically provided by 
reference views to (stable) cold space and onboard blackbody 
targets, whereas limb-scatter instruments either use the ex- 
traterrestrial solar spectra or limb measurements at upper tan- 
gent height(s) for normalization. Nadir viewing instruments 
provide full global coverage of the Earth on a daily basis. 


with 6 km vertical resolution for the stratospheric ozone pro- 
file, and between tens and hundreds of thousands of profiles 
per day. Detailed information on ozone data collected by 
limb sounders, their vertical and spatial resolution, and com- 
parisons between individual ozone climatologies for long- 
term changes is provided in a study by Tegtmeier et al. (2013) 
about the Stratospheric Processes and their Role in Climate 
(SPARC) initiative on ttace gas and aerosol climatologies. 

Temporal coverage of each measurement system is shown 
in Fig. 2. Average spatial coverage for the solar-occultation 
instruments is shown in Fig. 4 since this is very distinct 
for each of the systems. Spatial coverage for the stellar- 
occultation, limb and nadir systems described in the follow- 
ing sections is shown in detail in Toohey et al. (2013), and 
is therefore not explicitly repeated here. Figure 3b-d show in 
detail the vertical resolution of each of the satellite measure- 
ment systems that are described below. 

3.1 Solar-occultation instruments 

3.1.1 SAGE 

All of the Stratospheric Aerosol and Gas Experiment 
(SAGE) instruments (SAGE I on the Applications Explorer 
Mission-B (AEM-B), SAGE II on the Earth Radiation Bud- 
get Satellite (ERBS) and SAGE III on the Meteor-3M space- 
craft), used solar occultation to take spectral ttansmission 
measurements along a slant path through the atmosphere. 
SAGE III added lunar occultation measurements. The SAGE 
instruments scanned a small field of view across the solar 
disk every 2 s. The wavelength-dependent atmospheric trans- 
mission is determined by taking the ratio of a measurement 
looking through the atmosphere to one taken when the sun or 
moon is above the atmosphere as seen from the spacecraft. 
The systematic uncertainty of the Uansmission is indepen- 
dent of instrument characteristics as long as the instrument is 
stable over the brief few minutes spanned by each occultation 
event. Ozone is inferred from spectral measurements near 
600 nm, at the peak of the Chappuis band and nearby chan- 
nels are used primarily to characterize the spectral aerosol 
extinction. The small field of view and high signal-to-noise 
ratio enabled by looking at such bright sources enables the 
retrieval of the ozone profile at 1 km vertical resolution 
with low random uncertainty (~ 1 %) from just below the 
Uopopause to well above the stratopause (Table la, b). 

SAGE I made measurements over a 33 month period, start- 
ing in February 1979 (Fig. 2). SAGE II made measurements 
for 21 years starting in October 1984 (Table lc) covering the 
time when ozone depletion was at its peak as well as dur- 
ing the transition to lower chlorine loading in response to 
the Montreal Protocol. SAGE III made measurements from 
December 2001 to December 2005; it provided complemen- 
tary high-latitude coverage (Fig. 4). 

The solar occultation technique produces a very stable 
long-term data set. The ratio between scans across the solar 
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Latitude-Time Sampling at 30km by Solar Occultation Instruments 



Figure 4. Latitudinal coverage for the solar-occultation instruments 
(example years: 2003 and 2005). SAGE I is similar to SAGE II 
(blue triangles), POAM and ILAS instruments are similar to SAGE 
III (turquoise triangles). 


disk above the atmosphere and scans through the atmosphere 
to determine the spectral extinction vs. altitude are calcu- 
lated. The instrumental radiometric properties divide out in 
this ratio. Scanning across, rather than staring at, the sun pro- 
vides accurate knowledge of the altitude registration of each 
measurement. 

The most recent version of the SAGE data is version 7. 
This is the first version that will unify the algorithm, an- 
cillary data and spectroscopy across all three SAGE mis- 
sions (Damadeo et al., 2013). Data are archived at the Atmo- 
spheric Science Data Center, NASA Langley Research Cen- 
ter in Hampton, Virginia (see Table Id for URL). 

3.1.2 HALOE 

The Halogen Occultation Experiment (HALOE) was 
launched aboard the Upper Atmosphere Research Satellite 
(UARS) in 1991 and operated successfully for fourteen years 
until November 2005. HALOE made global observations 
verifying the effects of chlorine compounds on the chemical 
loss of ozone in the upper stratosphere. In addition to HC1 
and ozone, HALOE provided profile measurements of HiO, 
CH 4 , NO, NO 2 , HF, temperature and aerosol extinction us- 
ing the technique of solar occultation. Its measurements of 
HF, HC1, CH 4 , and NO were obtained using gas filter corre- 
lation radiometry. HALOE ozone transmission profiles were 
measured by the more traditional method of broadband filter 
radiometry (Russell III et al., 1993), using a broadband chan- 
nel from 9.2 to 10.4 micrometers. Overall measurement un- 
certainties range from 9-25 % in the lower stratosphere and 
from 9-20 % in the upper stratosphere and mesosphere. 

The HALOE retrieval algorithm employed a top-down, 
onion-peeling approach and iterated the calculated tangent- 
layer transmissions to achieve a match with the measured 
transmissions from each of its channels. Measured transmis- 
sions are calibrated according to their solar look values at 
top of scan. The one exception was the method for obtaining 
the temperature profile, which began with the assignment of 
the co-located National Center for Environmental Prediction 
(NCEP) temperature and pressure (or T(p) value) at 3 1 .5 km 
to the HALOE-measured CO 2 channel transmission profile 


at that same level. Then, a first-guess T ( p ) profile was iter- 
ated and proceeded upward by matching the calculated trans- 
mission value in each tangent layer with that observed in 
the HALOE CO 2 channel centered at 2.8 pm. The measured 
transmission profiles for all the species, including ozone, 
were then registered vs. the pressure profiles obtained in this 
way, prior to final retrievals of the species volume mixing ra- 
tio profiles (Table la-c). The HALOE Level 2 profiles can be 
downloaded from the NASA Goddard Earth Sciences Data 
& Information Services Center (GES DISC) or the HALOE 
webpage (see Table Id for the URLs). 

3.1.3 ILAS and ILAS-II 

The Improved Limb Atmospheric Spectrometer (ILAS) 
was a satellite-borne solar-occultation sensor on board the 
Japanese Advanced Earth Observing Satellite (ADEOS) 
(Sasano, 2002). ILAS consisted of an infrared spectrometer 
that covers the wavelength region from about 6 to 12 pm with 
a detector array of 44 elements and a visible spectrometer 
from 753 to 784 pm with a detector array of 1024 elements. 
ADEOS was successfully launched in a sun-synchronous or- 
bit on 17 August 1996. After a 3 month initial checkout pe- 
riod, continuous operation of ILAS started on 30 October, 
lasting until 29 June 1997 when a failure in the satellite’s so- 
lar battery system occurred (Fig. 2). For a period of about 8 
months, ILAS made solar-occultation measurements at 57- 
71° N and 64-88° S, and collecting ^6700 profiles of O 3 , 
HNO 3 , NO 2 , N 2 O, CH 4 , H 2 O and aerosol extinction coeffi- 
cients at 780 nm. 

ILAS-II was a satellite-borne solar-occultation sensor on 
board the satellite ADEOS-II (Nakajima, 2006). ILAS-II 
consisted of four spectrometers: an infrared spectrometer 
(between 6.21pm to 11.76 pm); a mid-infrared spectrome- 
ter (between 3.0 pm and 5.7 pm); a high-resolution infrared 
spectrometer (between 12.78 pm and 12.85 pm); and a visible 
spectrometer (between 753 nm and 784 nm). ADEOS-II was 
successfully launched on 14 December 2002, also in a sun- 
synchronous orbit. After a 3 month initial checkout period, 
continuous operation of ILAS-II started on 2 April 2003, 
lasting until 24 October 2003 when the satellite’s solar bat- 
tery system malfunctioned (Fig. 2, Table lc). For a period 
of about 7 months, ILAS-II made solar-occultation mea- 
surements at 57-71° N and 64-88° S, and gathered 5700 
profiles of O 3 , HNO 3 , NO 2 , N 2 0, CH 4 , H 2 0, CIONO 2 , 
N 2 O 5 , CFC-11, CFC-12 and aerosol extinction coefficients 
at 780 nm. 

Ozone and other minor gas vertical profiles were retrieved 
by applying both the onion-peeling method and the nonlin- 
ear least-squares fitting method. Version 5.20 ILAS and ver- 
sion 1.4 ILAS-II ozone profiles were validated by comparing 
with various data sources (Sugita et al., 2002, 2006). These 
validations showed that ILAS ozone data agreed with other 
data sets within ± 10 % with a few exceptions between 11 
and 64km (Table la, b). ILAS-II ozone data agreed with 
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other data sets within ±10% at altitudes between 11 and 
40 km for the Northern Hemisphere. Above 41 km, a nega- 
tive bias between the NH ILAS-II ozone data and the other 
data increased, reaching 30% at 61-65 km. Southern Hemi- 
sphere ILAS-II ozone data agreed with other data sets within 
±10% at altitudes between 11 and 60 km. Since the solar- 
occultation technique is used with both ILAS and ILAS-II, 
the instruments are self-calibrated by measuring 100% and 
0 % intensity of the sun each time before and after the atmo- 
spheric measurement. ILAS ozone data are archived at the 
ILAS and ILAS-II webpages (see Table Id for URLs). 

3.1.4 POAM II and POAM III 

The Polar Ozone and Aerosol Measurement II (POAM II) 
instrument was developed by the Naval Research Labora- 
tory (NRL) to measure the vertical distribution of atmo- 
spheric ozone, H 2 O, NO 2 , aerosol extinction and temper- 
ature using with the solar-occultation technique. POAM II 
measured solar extinction in nine narrow band channels, cov- 
ering the spectral range from approximately 350 to 1060 nm. 
POAM II was launched aboard the French SPOT-3 satellite 
on 26 September 1993 into a sun-synchronous polar orbit. 
The POAM II mission was ended by the failure of the SPOT- 
3 satellite in November of 1996 (Fig. 2). 

The POAM III instrument was a satellite-based nine- 
channel, visible-near-infrared photometer for making mea- 
surements of ozone, aerosol extinction, H 2 O and NO 2 in 
the polar stratosphere using the solar occultation technique 
(Lucke et al., 1999). The POAM III instrument was devel- 
oped by NRL and measured solar extinction in nine narrow 
band channels, covering the spectral range from approxi- 
mately 350 to 1020 nm. POAM III was carried by the French 
SPOT-4 spacecraft, which launched on an Ariane-4 rocket in 
March 1998. The POAM III instrument failed on 5 Decem- 
ber 2005. 

POAM made 14—15 measurements per day in each hemi- 
sphere around a circle of latitude. In the Northern Hemi- 
sphere the measurement latitude ranged from 55° N to 73° N, 
and in the Southern Hemisphere, 63° S to 88° S. The POAM 
III retrieval version 4 includes the following profile data: 
ozone concentration (5-60 km), NO 2 concentration (20- 
40 km), aerosol extinction at 6 wavelengths (5-25 km), H 2 O 
(5-45 km) and UK Met Office temperature, pressure, and po- 
tential vorticity (0-60 km), all with a 1 km vertical resolu- 
tion. POAM II has a slightly smaller altitude range for the 
retrieved profiles than POAM III (Table la, b). The estimate 
for random error on the retrieved concentrations is approxi- 
mately 3-5 % above 20 km and increasing at lower altitudes. 
For a complete error analysis see Lumpe et al. (2002). As so- 
lar occultation instruments, both POAM systems are largely 
self-calibrating. No drift of degradation of the POAM III 
ozone retrievals has been observed. 

The most recent version of the POAM retrieved profiles 
is version 4 (pdf document on the POAM ftp archive, see 


Table Id for URL), which differs only very little from version 
3 (Lumpe et al., 2002). Data is available either through the 
NASA Langley Atmospheric Science Data Center or through 
the Naval Research Laboratory (see Table Id for URLs). 

3.1.5 ACE-MAESTRO 

The atmospheric chemistry experiment (ACE) Measurement 
of Aerosol Extinction in the Stratosphere and Troposphere 
(MAESTRO) occultation instrument (McElroy et al., 2007; 
Nowlan et al., 2007), on board the Canadian SCISAT satel- 
lite, is a dual, diode-array spectrometer measuring in the UV- 
VIS-NIR spectral range with a nominal wavelength range of 
285-1015 nm (525 and 1010 nm for its primary ozone mea- 
surement mode in solar occultation). SCISAT was designed 
to measure the composition of the Earth’s atmosphere from 
space, with particular emphasis on the middle atmosphere 
ozone distribution, and related trace gases in the Arctic. The 
satellite was launched on 12 August 2003 and designed for 
a 2 year mission and is currently still in operation (Fig. 2). It 
is in a 74° inclined circular orbit at an altitude of 650 km. 

MAESTRO makes measurements in the solar occultation 
mode and shares a common sun tracker and optical view- 
ing direction with the atmospheric chemistry experiment- 
Fourier transform spectrometer instrument on the same satel- 
lite. Therefore, they observe almost exactly the same slant 
column of air, albeit with different fields of view, and with 
only a slight difference resulting from differential refractive 
effects. MAESTRO has a relatively small (0.02°) optical res- 
olution which provides a resolution of~ 1 km in altitude (Ta- 
ble lb). 

The primary purpose in adding MAESTRO to the satel- 
lite payload was to obtain wavelength-dependent, aerosol ex- 
tinction information. However, in order to remove the ef- 
fects of gaseous absorbers, a spectral fitting analysis is done 
which produces ozone, NO 2 and water vapor slant column 
amounts. Vertical profiles of these species in volume mix- 
ing ratio are also produced, covering the altitude range of 
5-35 km (Table la). Overall uncertainty on the retrieved pro- 
files is 10 %, with smaller uncertainties («a 3 %) in the mid- 
to upper stratosphere (Dupuy et al., 2009; Kar et al., 2007). 

Because MAESTRO is measuring in occultation mode, 
a reference spectrum is collected before or after each oc- 
cultation. Therefore, the optical depths calculated are almost 
independent of instrument sensitivity. The UV half of MAE- 
STRO has lost sensitivity steadily throughout the mission, 
with the loss of sensitivity migrating from shorter to longer 
wavelengths to the point where the instrument has produced 
no useful data for the last five years. However, the visible- 
near-IR spectrophotometer has shown very little change in 
sensitivity. Based on the behavior of the two instruments in 
the overlap region of the spectrum it is likely that the loss 
of sensitivity is due to attenuation taking place outside the 
instrument. 
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MAESTRO data can be obtained from the ACE Science 
Team at the University of Waterloo, Canada (see Table Id 
for URL). 

3.1.6 ACE-FTS 

The atmospheric chemistry experiment-Fourier transform 
spectrometer (ACE-FTS), on board the SCISAT satellite, is 
a solar occultation instrument measuring in the mid-infrared 
spectral region (Bernath, 2006). SCISAT began routine mea- 
surements on 2 1 February 2004 and is still currently in oper- 
ation (Fig. 2). 

ACE-FTS is a high-resolution (0.02 cm -1 ) FTS that mea- 
sures the full spectral range between 750 and 4400 cm -1 
(Bernath et al., 2005). Each day ACE-FTS measures approx- 
imately 15 sunrise and 15 sunset occultations and achieves 
nearly complete global latitude coverage over a period of 
three months (e.g., one “season”; Table lc). These mid- 
infrared spectral measurements extend from the cloud tops 
to 150km and ozone is retrieved up to 95km (Table la). 
The vertical spacing between each 2 s ACE-FTS measure- 
ment varies between 1.5 and 6 km depending on the satel- 
lite’s orbit geometry. The vertical sampling is 3-4 km based 
on the field of view of the instrument (Table lb). 

Constituent profiles are retrieved from the ACE-FTS trans- 
mission spectra as volume mixing ratios (VMRs) using 
an unconstrained non-linear least-squares global fitting ap- 
proach (Boone et al., 2005). Initially, CCU transitions are 
used to determine the pressure and temperature as a func- 
tion of altitude between 12 and 120 km. Below 12 km, mete- 
orological results from the Canadian Meteorological Centre 
(CMC) operational weather analysis and forecast system are 
used (note that in version 3.5 of the ACE-FTS retrievals, the 
CMC pressure and temperature results are used below 15 km, 
Boone et al., 2013). The resulting temperature and pressure 
profiles are the basis for retrieving VMR profiles of ozone 
and over 30 trace gas species. For the ozone retrieval, spec- 
tral features between 984 and 1 190 cm -1 are used. 

The spectroscopic parameters on which the calculations 
are based are derived from the High-resolution TRANsmis- 
sion (HITRAN) molecular absorption database 2004 line 
list (Rothman et al., 2005). The version 2.2 ozone updated 
(herein v2.2) data set was validated extensively by Dupuy 
et al. (2009). In comparisons with ozonesondes and satel- 
lite limb sounding instruments, a small positive bias of + 1 to 
+8% (+5% on average) was found between 16 and 44 km. 
At altitudes above 45 km, larger differences of up to +40 % 
(+20% on average) are found for the v2.2 ozone profiles. 
The current version of the ACE-FTS data set is version 3.5 
(Boone et al., 2013). Waymark, et al. (2013) have evalu- 
ated the differences between the well-validated v2.2 and new 
v3.5 profiles and found that there is a reduction of ~ 5 % 
at altitudes above 35 km in v3.5 compared to v2.2. There is 
no detailed error budget available for the ACE-FTS ozone 


product at present. However, there is an associated fitting un- 
certainty provided for each measurement. 

The ACE-FTS instrument is self-calibrating. A series 
of reference exo-atmospheric spectra (“high sun” measure- 
ments) are collected during each orbit and are used with 
the atmospheric measurements to calculate the transmittance 
spectra that are then used in the retrieval process. Data are 
made available by the ACE team on request or via web- 
page of the ACE Science Team at the University of Waterloo, 
Canada (see Table Id for more detail). 

3.2 Stellar occultation instruments 

3.2.1 GOMOS 

The Global Ozone Monitoring by Occultation of Stars (GO- 
MOS) instrument is one of the three atmospheric instru- 
ments on-board the European Space Agency’s Environmen- 
tal Satellite (Envisat) satellite. The satellite was launched in 
March 2002 into a polar sun-synchronous orbit with a de- 
scending node at 22:00 LT GOMOS measurements ceased 
in April 2012 when contact with Envisat was lost (Fig. 2). 
About 410000 nighttime ozone profiles were measured by 
GOMOS during its 10 year period of operation. 

The GOMOS spectrometer covers the wavelength range 
from ultraviolet, visible and near infrared (UV-VIS-NIR). 
It uses the stellar occultation technique for observing trace 
gas and aerosol profiles from the upper troposphere to the 
mesosphere (Bertaux et al., 2010). Approximately 160 stars 
are bright enough for atmospheric profiling, providing global 
coverage from pole to pole. The altitude range is typically 
from the top of the clouds to 100km (Table la). Geophysi- 
cal validations have shown that the nighttime measurements 
of ozone are accurate with a ±2 % bias in the stratosphere 
between 20-40 km (e.g. Meijer et al., 2004; van Gijsel et al., 
2010). The sampling resolution of GOMOS is 0.5-1 .7 km de- 
pending on the geometry of the occultation (Table lb). Due 
to refraction the sampling resolution improves below 40 km. 
The stars are point sources, which ensures very good knowl- 
edge of the measurement position (< 30 m). 

The retrieval of ozone and other trace gases using 
stellar occultation method is relatively straightforward 
(Kyrola et al., 2010). The present retrieval assumes known 
Rayleigh extinction, which is estimated using the European 
Centre for Medium-Range Weather Forecasts (ECMWF) and 
MSIS90 models for neutral density profiles. Only a mini- 
mal amount of a priori information on ozone is needed and 
it is included in the retrieval by using Tikhonov regulariza- 
tion method assuming smoothness of the ozone profiles. The 
vertical resolution of the ozone profiles is thus lower than the 
actual sampling resolution: 2 km in the stratosphere and 3 km 
in the mesosphere. Ozone is retrieved together with aerosols, 
NO2, and NO3 using wide wavelength band 250-675 nm. 
The stellar characteristics (temperature and brightness) are 
linked to signal-to-noise ratio resulting in a varying random 
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uncertainty of the data depending on the stars used. The ran- 
dom uncertainty of ozone measurements in the stratosphere 
ranges from 0.5-4 %, except in the lowermost stratosphere 
where values are about 10 %, and 2-10 % in the mesosphere. 
The largest systematic error contribution is due to uncertain- 
ties in aerosol modeling and is about 20% below 20 km, 
1-5 % at 20-25 km and negligible above 25 km (Tamminen 
et al., 2010). 

As an occultation instrument GOMOS is, in principle, less 
sensitive to instrument degradation. However, careful eval- 
uation of possible drifts is needed since the signal to noise 
ratio is low with dim stars and the increasing amount of 
dark charge is known to cause problems despite the self- 
calibrating measurement principle. These problems are mini- 
mized by careful screening of the data for dim and cool stars. 

In addition to the nighttime occultation measurements dis- 
cussed above, GOMOS has daytime occultation measure- 
ments (Table lc). Presently GOMOS daytime occultations 
are of poor quality (e.g. van Gijsel et al., 2010). Daytime 
ozone profiles can, however, be derived from GOMOS mea- 
surements of scattered solar light (Tukiainen et al., 2011). 
These data are processed by the Finnish Meteorological In- 
stitute (FMI) and ozone profiles are available from the ESA 
SPIN project. All occultation measurements are available via 
the ESA EARTHNET webpage (see Table Id for more de- 
tail). 

3.3 Limb emission instruments 
3.3.1 LIMS 

The Limb Infrared Monitor of the Stratosphere (LIMS), 
an infrared, limb-emission instrument, operated on Nimbus 
7 following its launch in late October 1978 through late 
May 1979, at which time its solid cryogen coolant was de- 
pleted (Fig. 2). The quality of the original LIMS Version 5 
data set is described in Gille and Russell III (1984) and refer- 
ences therein. The LIMS retrieval algorithm uses a top-down 
onion-peeling approach and iterates to achieve a match of the 
calculated and measured limb radiances for its wide and nar- 
row C02 radiometer channels in the 15-micrometer region. 
A constant CO 2 volume mixing ratio profile of 325 ppmv is 
assumed for the radiance calculations in the upper strato- 
sphere, where the forward models for the wide and narrow 
CO 2 channels are aligned and iterated. Radiance profiles for 
the LIMS species channels are registered against the pressure 
profiles that are obtained in this manner. 

In order for the LIMS data set to be more compatible 
with those from instruments of the UARS, Envisat, and the 
EOS satellite programs of more recent decades, LIMS pro- 
files were retrieved anew using updated spectral line param- 
eters for CO 2 and other species, including ozone, from the 
1992 or 1996 HITRAN line lists (Remsberg et al., 2004). 
LIMS ozone profiles were obtained from infrared limb ra- 
diance measurements at 9 to 10 micrometers. The resultant 


LIMS Version 6 Level 2 profiles of temperature, species and 
height are tabulated at 18 levels per decade of pressure or at 
a spacing of 0.88 km (Table la). They have an effective ver- 
tical resolution of 3.7 km, in order to account properly for 
the effects of the instrument on the observed radiances and 
of the retrieved temperatures on the retrieved ozone profiles 
(Table lb). The Version 6 data set was processed in 2002 
and the Level 2 and 3 products are archived at NASA God- 
dard (GES DISC; see Table Id for URL). A more detailed 
description of the LIMS V6 ozone and its quality is provided 
in Remsberg et al. (2007). Overall uncertainties for LIMS 
ozone range from 20-30 % in the lower stratosphere, to ap- 
proximately 12 % in the mid- to upper stratosphere. 

Each measured, up/down horizon radiance scan pair in- 
cluded a view of its warm, internal blackbody, followed by 
a calibration measurement of cold space. Characteristics of 
the blackbody remained stable throughout the 7.25 months 
of the LIMS mission. 

Most recently, Remsberg et al. (2013) reported on the as- 
similation of the LIMS V6 ozone into the GEOS -5 model 
of NASA’s Modern-Era Retrospective Analysis for Research 
and Applications (MERRA) system. They reported a prob- 
lem with the LIMS ozone data at pressures higher than about 
50 hPa (volume mixing ratios less than 0. 1 ppmv at low lati- 
tudes) and, therefore, recommend that these data be screened 
out prior to using the LIMS data set for detailed scientific 
studies of lower stratospheric ozone. 

3.3.2 UARS MLS 

The Upper Atmosphere Research Satellite (UARS) Mi- 
crowave Limb Sounder (MLS) was one of ten instruments 
on the NASA UARS platform, launched on 12 Septem- 
ber 1991 (Reber et al., 1993, Fig. 2). UARS-MLS (a pre- 
decessor to Aura MLS) pioneered microwave limb sounding 
of the Earth’s stratosphere and mesosphere from space (see 
Barath et al., 1993; Waters, 1993). UARS-MLS measured 
microwave emission near 205 GHz for its primary strato- 
spheric ozone retrievals. The UARS orbit was inclined at 
57° and MLS data coverage was from 34° on one side of 
the equator to 80° on the other side, switching sides after 
each satellite yaw maneuver (which occurred roughly every 
36 days). The measurements thus cover all local solar times 
during each “yaw period”. Profiles were spaced ~ 3^f° along 
the orbit track and the average daily sampling in longitude 
is 12°; there are about 1350 profiles per day (Table lc). 
The frequency of UARS-MLS operational days generally de- 
creased from about 100 % in late 1991 through 1993, to about 
50% in late 1994 and down to as low as a few tens of days 
per year from 1995 onward. 

The UARS MLS data are considered most robust for time 
series analyses through mid-June 1997, after which cessation 
of 63 GHz observations (to save spacecraft power) resulted 
in the loss of retrieved temperatures. The retrieval algo- 
rithms were based on optimal estimation, as for Aura MLS, 
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but without two-dimensional retrieval capabilities to account 
for horizontal volume mixing ratio gradients. Version 5 was 
the most recent public release of UARS MLS O3 data (see 
Livesey et al., 2003) and the work by Froidevaux et al. (1996) 
on an earlier data version for more details and validation re- 
sults). 

UARS MLS ozone profiles are recommended from 100 to 
0.22 hPa and are retrieved on the same vertical grid as for 
Aura MLS v2.2 O3; the vertical resolution is 3.5-5 km in the 
stratosphere and 5-8 km in the mesosphere (Table la, b). Es- 
timated random errors are 0. 2-0.3 ppmv and typical system- 
atic error estimates (2er) are ~ 6 % , degrading to 15% (or 
more) near 68-100 hPa, where some high biases vs. correla- 
tive data were noted. The official public distribution location 
for UARS MLS data is at the NASA GES-DISC Mirador 
website (see Table Id for URL). 

3.3.3 Aura MLS 

NASA’s Earth Observing System (EOS) Microwave Limb 
Sounder (MLS) instrument was launched on the Aura satel- 
lite on 15 July 2004 (Fig. 2). This improved sounder (and 
successor to UARS MLS) measures a variety of atmospheric 
species (and temperature); O3 retrievals are provided from 
the upper troposphere to the mesosphere. Aura MLS per- 
forms limb scans ahead of the Aura satellite, which is in 
a sun-synchronous, near-polar orbit. Aura MLS measures 
ozone microwave emission for its primary (standard) prod- 
uct near 240 GHz. Retrieved profiles occur at the same lat- 
itude every orbit and are separated by & 170 km along the 
sub-orbital track. This yields nearly 3500 profiles every day 
(day and night) from 82° S to 82° N (Table lc). Waters et al. 
(2006) describe the instrument and its measurements. 

The Aura MLS retrievals are obtained via optimal esti- 
mation (Rodgers, 1976, 2000). Gauss-Newton iteration is 
used, with a second-order Tikhonov constraint (Tikhonov, 
1963). This constraint is applied to the profile second 
derivatives (vertically and horizontally). Specific retrieval as- 
pects include adaptation to a two-dimensional system, us- 
ing the line-of-sight measurements from several scans to 
derive multi-profile information (Livesey and Read, 2000). 
The Level 2 data include screening flags regarding instru- 
ment/retrieval status, quality of fit and retrieval convergence 
(see Froidevaux et al., 2008). Retrieval details are provided 
by Livesey et al. (2006) and references therein. 

Version 2.2 O3 is retrieved on a pressure grid with 6 lev- 
els per decade up to 0.1 hPa (and coarser vertical resolution 
above this); v3.3 data are provided on a finer grid and down to 
261 hPa, but are not recommended overall over v2.2 because 
of larger vertical oscillations in the upper troposphere and 
lower stratosphere (UTLS) at low latitudes (Table la). Aura 
MLS v2.2 single profiles are recommended for scientific use 
in the 215 hPa to 0.02 hPa range; the vertical resolution is 
ss 3 km (or slightly less) in the stratosphere, degrading to 4— 
6 km for pressures of 0.1 hPa or less (Table lb). Estimated 


random errors are 0.1 to 0.3 ppmv («a 2-15 %) and systematic 
error estimates (2er) are 5-10 % for most of the stratosphere 
and lower mesosphere (for more details and validation, see 
Froidevaux et al., 2008, and references therein). The offi- 
cial public distribution location for Aura MLS data is at the 
NASA GES-DISC Mirador website (see Table Id for URL). 

3.3.4 MIPAS 

The Michelson Interferometer for Passive Atmospheric 
Sounding (MIPAS) is a limb emission, Fourier transform 
spectrometer (Fischer et al., 2008). Looking horizontally 
through the atmosphere at various elevation angles, it mea- 
sured high-resolution, infrared emission spectra of atmo- 
spheric constituents, thus providing altitude -resolved infor- 
mation on these constituents. MIPAS was one of the atmo- 
spheric chemistry instruments on the ESA’s Envisat research 
satellite, which was launched into its sun-synchronous polar 
orbit on 1 March 2002. In March 2004, problems with the 
interferometer slide occurred and operation was interrupted; 
however, on 27 January 2005 operation in the nominal mode 
was resumed with reduced spectral but improved spatial res- 
olution (optimized resolution measurements). MIPAS pro- 
vided data until 8 April 2012 when contact to the satellite 
was lost (Fig. 2). 

There exist multiple MIPAS level-2 retrievals. One version 
was generated with the research processor developed and op- 
erated by the Institute of Meteorology and Climate Research 
(IMK), KIT, Karlsruhe, Germany, in cooperation with the In- 
stitute de Astrofisica de Andalucfa (IAA), CSIC, Granada, 
Spain (von Clarmann et al., 2003, 2009). The retrieval was 
performed on a fixed altitude grid usually finer than the tan- 
gent altitude spacing, and was regularized by a Tikhonov first 
order approach that leads to smoothing of the profiles without 
bias towards an a priori. IMK high-spectral-resolution ozone 
data (2002-2004) were validated by Steck et al. (2007) and 
the first validation of the reduced spectral resolution data was 
published by Stiller et al. (2012). 

Frequent gain and offset calibration by means of deep 
space and blackbody measurements minimize drifts of MI- 
PAS measurements. There is a higher order effect due 
to slowly decreasing detector non-linearity that is over- 
estimated towards the end of the MIPAS measurement 
period. This explains a predominantly negative drift of 
MIPAS-IMK ozone vs. MLS-Aura ozone, which can reach 
—0.2 ppmv decade -1 around 40 km (Eckert et al., 2014). 

Another MIPAS processor is run by ESA (Ridolfi et al., 
2000; Raspollini et al., 2006, 2013b). It is based on the algo- 
rithm ORM (optimized retrieval model) developed by an Eu- 
ropean Consortium led by IFAC-CNR (Istituto di Fisica Ap- 
plicata “Nello Carrara”, part of the Consiglio Nazionale delle 
Ricerche). This processor was specifically designed for oper- 
ating in near-real time, and hence for working automatically 
in different atmospheric conditions, and uses the minimum 
amount of a priori information that may introduce a bias 
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in the profiles. To this purpose, it uses the Levenberg- 
Marquardt regularization inside the iterations and an a pos- 
teriori Tikhonov regularization with a self-adapting con- 
straint that ensure that minimum constraints are used and the 
best vertical resolution obtainable from the measurements is 
achieved in all atmospheric conditions. ESA high-spectral- 
resolution ozone data (2002-2004) were validated by Cortesi 
et al. (2007). 

A third MIPAS retrieval algorithm was developed at 
Oxford University (Moore and Remedios, 2010) and is 
based on the optimal estimation approach (Rodgers, 2000). 
A fourth retrieval algorithm, GMTR, has been developed as 
a joint effort of the University of Bologna and ISAC-CNR 
(Carlotti et al., 2006). The algorithm performs simultaneous 
retrievals of pressure, temperature, water and ozone analyz- 
ing all the measurements acquired during a full Envisat or- 
bit using a two-dimensional approach. The GMTR algorithm 
has been used to analyze the full MIPAS mission and the 
results are stored into the MIPAS2D database (Dinelli et al., 
2010, 2012). All these algorithms use the same level- lb data, 
which is provided by ESA, but differ in some assumptions for 
the forward model and for the spectral intervals used in the 
analysis (Laeng et al., 2014; Raspollini et al., 2013a). 

Recently the full validation of the ozone products of the 
four cited algorithms (Laeng et al., 2014), an assessment of 
the long-term stability, and an assessment of the 10 year trend 
of the IMK/1AA data product (Eckert et al., 2014) have been 
completed. Overall uncertainties on the ozone profiles range 
from 35-40 % at 10 km to 8-1 1 % at around 30 km. Vertical 
resolution for optimized resolution phase (years 2005-2012) 
ranges from 2-2.4km at 10km and 3.2-3. 6km at 30km to 
5 km at 70km. Vertical resolution for full resolution phase 
(years 2002-2004) ranges from 4 km around 10 km to 5 km 
around 30 km, to 8 km at 70 km (Table lb). In order to avoid 
artifacts in trend estimation from two data subsets of differ- 
ent characteristics resulting from the different measurement 
modes before and after 2004, a dedicated scheme has been 
developed that is immune to biases between data subsets (von 
Clarmann et al., 2010). The IMK/IAA MIPAS data are ac- 
cessible at the IMK webpage (see Table Id for URL). The 
ESA MIPAS data are accessible at the ESA webpage (see 
Table Id for URL). The MIPAS2D database is undergoing 
a validation process (Casadio et al., 2013), but the data is ac- 
cessible through the ISAC-CNR web site (see Table Id for 
URL). The MIPAS data retrieved by the Oxford group is not 
publicly available at this point. 

3.3.5 Odin-SMR 

The Odin satellite orbit is approximately 600 km and sun- 
synchronous. The orbit has precessed over the eleven years 
since Odin was launched but the local times of the mea- 
surements have remained relatively stable (equator cross- 
ing times between 06:00-07:00 and 18:00-19:00). The Odin 
sun-synchronous orbit leads to daytime observations in the 


summer hemisphere and nighttime observations in the win- 
ter hemisphere. 

The Sub-Millimeter Radiometer (SMR) on board the Odin 
satellite, launched in 2001, makes limb measurements of 
thermal emission lines of ozone in several different spectral 
bands. The 501.8 GHz version 2.1 level-2 product provides 
stratospheric ozone data in the 12-60 km range (above 17- 
18 km at mid-latitudes) with 2.5-3. 5 km vertical resolution 
and single-profile random uncertainty of about 20 % (Jones 
et al., 2007; Jegou et al., 2008, Table la, b). The systematic 
error is estimated to be smaller than 0.75 ppmv (Urban et al., 
2005, 2006). 

Measurements in this observation mode were carried out 
on every 3rd day until April 2007 and on every other day 
thereafter (Fig. 2). The thermal emission technique allows 
ozone to be measured during day and night and global fields 
between ~ 83° S and ~ 83° N are produced during typically 
14—15 orbits per observation day based on up to 65 limb- 
scans per orbit. 

The calibration for SMR is done with a standard hot-cold 
radiance calibration scheme, based on the internal hot load 
(viewed at the end of each scan, i.e. every 90s) and cold sky 
(viewed each tangent-view, i.e. every 0.85 s). This leads to an 
excellent long-term stability of the measurements. There are 
no known signs of instrument degradation so far. 

Individual ozone profiles of the 501.8 GHz band are quite 
noisy since they are derived from a relatively weak line, but 
averages agree well with correlative measurements (Jones 
et al., 2009). Other SMR ozone products such as from the 
strong ozone line in the 544.6 GHz band (version 2.0) are 
less noisy but have larger (known) biases. Data are available 
from the Odin-SMR webpage (see Table Id for URL). 

3.3.6 SMILES 

The Superconducting Submillimeter-Wave Limb-Emission 
Sounder (SMILES) was developed and operated under the 
cooperation of the Japan Aerospace Exploration Agency 
(JAXA) and the National Institute of Information and Com- 
munications Technology (NICT) (Kikuchi et al., 2010). 
SMILES conducted its atmospheric observations from 
12 October 2009 to 21 April 2010 (Fig. 2) on the Japanese 
Experiment Module (JEM) aboard the International Space 
Station (ISS). SMILES used three observation bands in 
the frequency regions of 624.32-626.32 GHz and 649.12- 
650.32 GHz. This instrument measured about 1600 profiles 
per day (both daytime and nighttime) within the latitude 
range of 38° S-65° N. Owing to the non-sun-synchronous or- 
bit of the ISS, measurements of the atmosphere are taken at 
different local times. This data set provides unprecedented 
opportunities to study diurnal ozone variations throughout 
the stratosphere (Imai et al., 2013; Sakazaki et al., 2013; 
Kasai et al., 2013; Kreyling et al., 2013). 

The SMILES mission has two products for ozone. One is 
the JAXA product (Imai et al., 2013) and the other is the 
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NICT product (Kasai et al., 2013). Version 3.0 (1 18-12-0602) 
of O 3 profiles in SMILES JAXA product are retrieved based 
on the combination of Optimal Estimation method (Rodgers, 
2000) with Tikhonov Regularization method (L2 constraint). 
Degradation of the instrument is corrected by an engineering 
process that uses some parameters obtained by on-orbit cal- 
ibrations to obtain brightness temperature spectra from raw 
data (Ochiai et al., 2013). The SMILES ozone product pro- 
vides ozone concentration as a value of volume mixing ratio. 
The nominal retrieved altitude range is from 8 to 120 km, 
and the typical “useful altitude range” is from 16 to 95 km 
(Table la). The retrieval vertical grid step is 2.5 km through- 
out the retrieval altitude range (Table lb). The information 
on the convergence status of iteration in retrieval process, the 
validity of the observation altitude range, and the FOV inter- 
ference onboard the ISS, are represented as the “screening 
flags” and these flags are stored in the “status field” of each 
HDF file. 

The NICT products of ozone is called SMILES NICT 
Level-2 (L2) version 2.1.5 product, which used the version 
007 calibrated Level- lb (Lib) spectra. The detail of the eval- 
uation of the Lib spectrum, retrieval procedure including the 
theoretical error analysis, and comparison/evaluation of the 
ozone vertical profile was shown in Kasai et al. (2013). 

The public SMILES L2 and L3 (SMILES climatology) 
data are available from the Data ARchives and Transmission 
System (DARTS) in the Institute of Space and Astronautical 
Science (ISAS) of JAXA and from the NICT SMILES page 
(see Table Id for URLs). 

3.4 Limb scattering instruments 

3.4.1 SCIAMACHY 

The Scanning Imaging Absorption spectroMeter for Atmo- 
spheric CHartographY (SCIAMACHY) was a payload on 
Envisat launched in March 2002 (Burrows et al., 1995; 
Bovensmann et al., 1999, 2011). Due to the failure of com- 
munication to the Envisat platform, SCIAMACHY measure- 
ments ended in April 2012. SCIAMACHY is one of the new- 
generation space-borne instruments capable of performing 
spectrally resolved measurements in several different modes: 
alternate nadir and limb observations of the solar radiation 
scattered in the atmosphere or reflected by the Earth’s surface 
as well as observations of the light transmitted through the 
atmosphere during solar or lunar occultation. SCIAMACHY 
was a passive imaging spectrometer comprising eight spec- 
tral channels covering a wide spectral range from 214 to 
2386 nm. For the limb-scatter ozone profile retrieval selected 
wavelengths in the Hartley (ultraviolet) and Chappuis (visi- 
ble) ozone bands are used (Sonkaew et al., 2009). 

In the limb-scatter geometry, SCIAMACHY observed the 
atmosphere tangentially to the Earth’s surface starting at 

3 km below the horizon, i.e., when the Earth’s surface is 
still within the field of view of the instrument, and then 


scanned vertically up to the top of the neutral atmosphere 
(about 100km tangent height; Table la). The useful verti- 
cal range for ozone retrieval is about 10 to 60 km altitude 
and the vertical resolution of the retrieved ozone profiles is 
about 4 km (Table lb). In the nominal mode about 100 mea- 
surements per orbit with 14 orbits per day were performed. 
Global coverage of the sunlit part of the earth was achieved 
in 6 days (Table lc). 

Similar to other limb-scattering instruments, pointing un- 
certainty is a major error source (von Savigny et al., 2005). 
The systematic uncertainty of the pointing for the whole 
limb scan is estimated to be about 200 m. The relative 
pointing error between different tangent heights is negli- 
gible. Above 20 km measurements have uncertainties up 
to 15% (for a single profile), with a possible 8 % sys- 
tematic error. Extensive validation of the SCIAMACHY 
limb-scatter ozone profiles by comparisons with other satel- 
lites has been reported for an earlier retrieval version in 
Mieruch et al. (2012). A detailed error budget is provided in 
Rahpoe et al. (2013) and ozone trends during the last decade 
are reported in Gebhardt et al. (2014). Instrument calibration 
effects are largely removed by using ratios of the limb spec- 
tra at different tangent heights. SCIAMACHY version 2.5 
data are available from the webpage of the Institute of En- 
vironmental Physics, University of Bremen, Germany (see 
Table Id for URL). 

An alternative retrieval of ozone profiles from SCIA- 
MACHY limb observations is provided by the European 
Space Agency/DLR (Doicu et al., 2007). 

3.4.2 Odin-OSIRIS 

The Canadian-built, Optical Spectrograph and InfraRed 
Imaging System (OSIRIS) (Llewellyn et al., 2004) has 
been in routine operation on the Swedish spacecraft Odin 
(Murtagh et al., 2002) since the fall of 2001 (Fig. 2). OSIRIS 
measures vertically resolved profiles of spectrally dispersed 
limb scattered sunlight from 280 nm to 800 nm with approxi- 
mately 1 nm resolution. Typically the Odin spacecraft rotates 
to scan the OSIRIS single line of sight through a tangent al- 
titude range that extends from 10km to 60km (Table la). In 
the sunlit part of the orbit, a single vertical ozone profile is 
retrieved from a single Odin scan that typically takes 90s. 
This results in as many as 60 scans per orbit (Table lc). The 
Odin scan speed along with the OSIRIS imaging rate and in- 
stantaneous field of view combine to give an approximately 
1 .5 km sampling resolution in the UTLS and 2 km at higher 
altitudes (Table lb). 

The retrieval scheme used to infer vertical ozone profiles 
from limb-scattered sunlight was developed at the Univer- 
sity of Saskatchewan (Degenstein et al., 2009) and involves 
an adapted Multiplicative Algebraic Reconstruction Tech- 
nique and the SASKTRAN (Bourassa et al., 2007) radiative 
transfer model. The vertical ozone profile within the SASK- 
TRAN model is relaxed in a multiplicative fashion until the 
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model and measurements agree to within a certain tolerance. 
This method is extremely stable, and because the OSIRIS 
noise values are very low, the retrieved data quality is very 
good. Overall uncertainties range from 5 % between 25 km 
and 50 km, 15% below 20 km and 7 % between 50 km and 
55 km. 

Studies have shown that OSIRIS ozone compares well 
with coincident SAGE II ozone profiles with biases less than 
3 % from 18 to 50 km (Degenstein et al., 2009). OSIRIS re- 
trievals are self-calibrating as the radiance data is normalized 
against both higher altitudes and non-absorbing wavelengths. 
These normalizations stabilize the retrieval to the majority of 
known factors; the exception being upwelling radiation from 
the ground and lower troposphere. Data are available from 
the Odin-OSIRIS webpage (see Table Id for URL). 

3.5 Nadir-viewing instruments 

3.5.1 SBUV 

The Solar Backscatter Ultraviolet (SBUV) instruments are 
nadir-viewing instruments that infer total column ozone and 
the ozone vertical profile by measuring sunlight scattered 
from the atmosphere in the middle ultraviolet. The instru- 
ments use Ebert-Fastie type double monochromators and 
measure at twelve discrete wavelengths from 252-340 nrn. 
The ratios of the backscattered radiance to direct solar irra- 
diance are used to infer ozone concentration as a function 
of altitude. Optimal estimation techniques are used to invert 
the measurements and estimate ozone amounts as a function 
of pressure (Table la). The vertical resolution of SBUV is 
6 km in the upper stratosphere (~ 3 hPa) and gets progres- 
sively worse at higher and lower altitudes, particularly in the 
tropics (Table lb). Figure 9 in Bhartia et al. (2013) provides 
more details. 

The SBUV record is the longest available for satellite- 
based profile ozone measurements from a single instrument 
type (McPeters et al., 2013). The record began with the 
launch of Nimbus-4 BUV in 1970. An improved version, 
SBUV, was launched in October 1978 on the Nimbus-7 satel- 
lite. In 1984 the first of seven SBUV/2 instruments was 
launched on the NOAA-9 satellite, with subsequent launches 
on the NOAA 11, 14, 16, 17, 18, and 19 satellites (Fig. 2). 

The instruments are flown in sun-synchronous polar orbits 
to obtain global coverage. Since the measurements rely on 
backscattered sunlight, data are taken only on the daylight 
side of each orbit and no data are taken during polar win- 
ter. The NOAA polar orbiting satellites were not exactly sun- 
synchronous and consequently drifted from near-noon orbits 
with NOAAs 9-14 drifting the most. NOAA-1 1 for example, 
drifted from an initial 13:30 equator-crossing time to a 17:00 
orbit in six years, returning as a morning orbit for the later 
part of the data record. Starting with NOAA-16, the orbital 
drift was partially mitigated, providing ~ 5 years of data be- 
fore significant drift began. As the orbit drifts the solar zenith 


angle at which the measurements are made increases and lat- 
itude coverage decreases, which can be a problem for trend 
determination (Table lc). 

The most recent retrieval algorithm is Version 8.6 (Bhartia 
et al., 2013). Significant changes from Version 8 include new 
ozone (McPeters and Labow, 2012) and cloud climatologies 
(Joiner and Vasilkov, 2006), a change in ozone cross sections 
to Brion, Daumont, and Malicet, and updated instrument cal- 
ibration (DeLand et al., 2012). The absolute calibrations of 
the individual instruments are adjusted at the radiance level 
to improve the consistency over the multi-instrument pe- 
riod. This is accomplished using coincident measurements 
during periods of overlap referred to as the“no local time” 
correction. Local solar time coincidence is required as well 
as spatial coincidence to account for real diurnal variations 
in ozone (DeLand et al., 2012). The instruments have been 
cross-calibrated during two time periods. The calibration of 
NOAA-17 SBUV/2 is used as the standard for NOAA-14, 
NOAA-16, and NOAA-18, while the NOAA-1 1 base cali- 
bration is used as the reference for the Nimbus-7 and NOAA 
9 absolute calibrations. No instrument overlaps Nimbus-4 
BUV; its absolute calibration is based on comparisons with 
ground-based instruments and has not changed in Version 
8.6. 

The primary data product for time-series analysis is the 
monthly zonal-mean data, available at the NASA God- 
dard DISC webpage (see Table Id for URL). Validation 
of SBUV monthly zonal-mean profiles (Kramarova et al., 
2013b) shows that the mean biases and standard deviations 
are mostly within 5 % in the stratosphere between 25 and 

1 hPa. The SBUV vertical resolution decreases in the lower 
stratosphere and troposphere, and as a result the SBUV 
smoothing error, defined as the component of vertical ozone 
variability that the observation system cannot measure, in- 
creases from < 1 % in the middle stratosphere up to 10-15 % 
in the lower stratosphere and troposphere. Combining layers 
in the troposphere/lower stratosphere and above 1 hPa is rec- 
ommended to limit smoothing errors, as described in Kra- 
marova et al. (2013a). 

In order to make the SBUV measurements accessible for 
trend analysis, it is necessary to combine the data from the 
different NOAA satellites. Two versions of a “merged” data 
set are available (NASA Goddard, S. Frith, personal commu- 
nication, 2013; NOAA CPC, J. Wild, personal communica- 
tion, 2013). Both are described in more detail in the second 
and third SI 2 N overview papers (Lambert et al., 2014; Harris 
et al., 2014). 

3.5.2 GOME/GOME-2 

The Global Ozone Monitoring Experiment (GOME) is 
a nadir-scanning ultraviolet and visible (UV-VIS) spectrom- 
eter. Its launch on board the European Remote Sensing 

2 satellite in April 1995 started a new era of measuring 
Earth-reflected solar radiation with relatively high spectral 


Atmos. Meas. Tech., 7, 1395-1427, 2014 


www.atmos-meas-tech.net/7/1395/2014/ 



B. Hassler et al.: Ozone profile measurements 


1415 


resolution for UV-VIS wavelengths (Fig. 2). SCIAMACHY, 
GOME-2 and OMI followed the GOME heritage. The UV- 
VIS spectrum between 270 nm and 340 nm is used for de- 
riving global, height-resolved information on the ozone dis- 
tribution in the atmosphere (Munro et al., 1998). As for 
SBUV, Optimal Estimation techniques (Rodgers, 2000) are 
used for the retrievals. GOME ozone profiles were exten- 
sively analyzed by Meijer et al. (2006) and when the OMI 
measurements became available these were validated by 
Kroon et al. (201 1). 

A common feature of nadir UV-VIS ozone profiles is their 
limited vertical resolution, estimated typically at 7-15 km, 
with a considerable amount of vertical smoothing (Table lb). 
Smoothing of the retrieved ozone profile is described by the 
averaging kernel, which is part of the retrieval results. Un- 
certainty of the retrieved ozone profile is determined by mea- 
surement error, a priori error and by smoothing error (about 
30 % in the troposphere, and ss 10 % in the stratosphere). The 
ozone profile retrieval depends on the absolute calibration of 
the instruments and therefore it is sensitive to degradation ef- 
fects. An offline degradation correction is needed in case of 
degradation of the instrument. 

Within the climate change initiative (CCI) program of 
ESA, a new algorithm has been developed to generate high- 
quality data for nadir ozone profiles from GOME, SCIA- 
MACHY and GOME-2. Data are available from the EU- 
METSAT archive, and it is planned to make it accessible 
via CCI web portal (see Table Id for URLs). Furthermore, 
GOME-2 profiles are retrieved on an operational basis by 
KNMI within the framework of EUMETSAT’s Ozone and 
Atmospheric Chemistry Monitoring SAF (03MSAF) (e.g. 
Mijling et al., 2010; van Peet et al., 2009) that can be ac- 
cessed via the Eumetcast dissemination (see Table Id for 
URL). 

3.5.3 OMI 

The Ozone Monitoring Instrument (OMI) is an imaging 
nadir viewing UV-VIS spectrometer on board of the Aura 
satellite mission (Levelt et al., 2006a, b). Due to the ap- 
plication of two-dimensional detectors OMI combines daily 
global coverage with a relative high spatial resolution of 
13 km x 24 km at nadir for wavelengths from 308 to 500 nm 
and 13 km x 48 km for wavelengths below 308 nm. The OMI 
science data record currently extents from October 2004 
to present (Fig. 2). The OMI instrument is extremely sta- 
ble radiometric ally, with only a few percent of degrada- 
tion even in the shortest ultraviolet channels, making it 
a unique data set for ozone monitoring. Since the be- 
ginning of 2009 part the OMI swath is affected by the 
so-called row anomaly (http://www.knmi.nl/omi/research/ 
product/rowanomaly-background.php). 


The OMI instrument will be succeeded by the TROPOMI 
instrument on board of the Sentinel-5 Precursor satel- 
lite, planned for launch late 2015 (Veefkind et al., 2012). 
TROPOMI uses the same observation technique as OMI, and 
includes several improvements, including the signal-to-noise 
ratio, the spatial resolution and additional spectral bands. 
TROPOMI will have a separate spectral band from 270 to 
320 nm, which is expected to lead to improved ozone profile 
retrievals. 

The ozone profile retrieval for OMI is based on the opti- 
mal estimation method (Rodgers, 2000). As discussed above 
in the GOME/GOME-2 section (Sect. 3.5.2), the vertical in- 
formation is limited, especially in the UTLS and the tro- 
posphere (Table la-c). The OMI ozone profile data are 
available at the NASA Goddard DISC (see Table Id for 
URL). The validation status of this product is described in 
Kroon et al. (2011) based on comparisons with ozoneson- 
des and several satellite data products. Since the begin- 
ning of 2009 part of the OMI swath is affected by the 
so-called row anomaly (http://www.knmi.nl/omi/research/ 
product/rowanomaly-background.php). This has several ef- 
fects on the data quality and some of these depend on the po- 
sition in the orbit. Currently the effect of the row anomaly on 
the ozone profile quality in the Northern Hemisphere has not 
yet been fully quantified. In the future, further improvements 
to the OMI ozone profile product are planned using a new 
version of the retrieval algorithm including accurate flagging 
of the row anomaly. In addition to the product available at 
the NASA DISC, there is also a research product (Liu et al., 
2010), with a strong focus on the retrieval of tropospheric 
ozone information. 

3.5.4 IASI 

The Infrared Atmospheric Sounding Interferometer (IASI) 
instrument is a nadir-viewing, Fourier transform spectrom- 
eter. It was first launched on board the sun-synchronous po- 
lar orbiting MetOp-A satellite on 19 October 2006 (Fig. 2). 
A second instrument was launched on 17 September 2012 on 
board Metop-B and, together with the forthcoming Metop-C, 
will ensure products continuity for at least 15 years. 

IASI scans across track in 30 successive steps (four pix- 
els are acquired simultaneously) either side of nadir, with 
a total swath of around 2000 km (Table lc). This allows for 
global coverage twice a day. Each field of view has a typ- 
ical footprint size of 12 km at nadir. IASI measures spec- 
tra in the thermal infrared between 645 and 2760 cm -1 at 
a spectral sampling of 0.25 cm -1 (0.5 cm -1 resolution after 
apodization). The 8461 channels are distributed in 3 bands, 
which, due to their spectral coverage, allow for the retrieval 
of surface and cloud properties, temperature and humidity 
profiles as well as ozone and a number of other atmospheric 
constituents (e.g. CO, CH4, NH3, SO2). The ozone profiles 
are retrieved in near-real time, using the FORLI-03 (Fast 
Optimal Retrievals on Layers for IASI) software (Hurtmans 
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et al., 2012), and will be available as part of the IASI Level 2 
product processing facility (August et al., 2012, Table la-c). 
Validation was performed using ground-based (Anton et al., 
2011; Dufour et al., 2012), aircraft (Parrington et al., 2012; 
Pommier et al., 2012), balloon (Gazeaux et al., 2013) and 
other satellite data (Scannell et al., 2012). 

IASI is auto-calibrated, with high radiometric accuracy. 
Retrieval methods for profile measurements use standard 
Optimal Estimation Method (OEM) or regularization tech- 
niques. Overall uncertainties in the ozone profiles range from 
about 20-30 % in the altitude range 0-6 km, and are about 
10-20% in the stratosphere and upper troposphere. At the 
moment, IASI ozone data are available from the IASI team 
on request. Data should be available via Eumetcast in 2014 
(see Table Id for more detail). 

3.6 Additional instruments 

Beyond the instruments/missions described in detail above, 
several other satellite instruments also provide ozone data. 
These will only be mentioned briefly in the following para- 
graphs since these instruments were either not used in any 
SI 2 N activity (e.g. included in merged data sets, or updated 
and improved retrievals) or have their main focus on tropo- 
spheric rather than stratospheric measurements. However, for 
the sake of completeness we include them here. 

The TES (Tropospheric Emission Sounder) instrument on 
board of Aura is a nadir-looking high-resolution infrared- 
imaging Fourier transform spectrometer observing in nadir, 
which provides profiles of ozone and many other trace 
species in the troposphere and lowermost stratosphere. TES 
provides about 1000 observations a day, with a nadir foot- 
print of 5 km x 8 km. The sensitivity of the measurements 
is highest in the mid-troposphere. The data record has been 
thoroughly validated (Worden et al., 2007). 

The HIRDLS (High-Resolution Dynamics Limb Sounder) 
instrument is a 21 channel (from 6.12 to 17.76 microns) limb 
scanning infrared radiometer measuring emission from the 
Earth’s limb onboard Aura. These measurements are inverted 
to yield global distributions of temperature, clouds, aerosols, 
and concentrations of ozone and other trace gases at high 
vertical and horizontal resolution from the upper troposphere 
into the mesosphere. After launch, activation of the HIRDLS 
instrument revealed that the optical path was blocked so that 
20 % of the aperture could view the Earth’s atmosphere. Nev- 
ertheless, measurements at high vertical resolution could be 
made at one scan angle. Geographically the data extend from 
63° S to 80° N with 1° spacing, and cover the period from 
29 January 2005 to 17 March 2008, with very few gaps. The 
estimated total accuracy is 1-10%, varying with altitude and 
latitude, and 10 to > 50 % in the tropical UTLS (Nardi et al., 
2008). 

The Sounding of the Atmosphere using Broadband Emis- 
sion Radiometry (SABER) instrument provides volume mix- 
ing ratio profiles of ozone from 15-100 km retrieved from 


mid-infrared limb emission measurements (Russell III et al., 
1999). The validation of the vl.07 ozone profiles in the mid- 
dle and upper atmosphere (~ 15-70 km) are described in 
Rong et al. (2009). 


4 Summary 

This paper describes long-term ozone profile data derived 
from five of the most common ground-based systems and 19 
satellite-based instruments. The physical principles of mea- 
surements, specifics of data retrievals and uncertainties of 
measurements have been described in more detail in vari- 
ous publications and reports with relevant references given in 
this paper. Therefore, the goal of this paper is to provide an 
overview of well-established ozone records, while also offer- 
ing updated information on recently reprocessed and newly 
available data sources. The summary includes links to data 
archives, lists uncertainties of measurements and retrievals, 
provides spatial and temporal characteristics of the data sets 
described above and indicates the latest released data version 
(see Table Id, the numbers given are indicative or typical val- 
ues). This overview does not represent a complete list of all 
available ozone profile measurements, but aims to provide 
the community with information about the latest improve- 
ments and recalculations of the most commonly used ozone 
profile time series data that are most relevant in the context 
of trend analysis, bundled in one place for easy access and 
comparison. 

Although according to Figs. 2 and 3 (also Fig. 1 in Toohey 
et al., 2013), it appears that the ozone profile measure- 
ments are available globally for the entire 1960-2000 pe- 
riod, one needs to keep in mind that there are regions of 
the globe and times when data coverage is sparse. A small 
number of quasi-continuous ground-based measurements by 
the Umkehr method began in the early 1960s and only a few 
ozonesonde locations provide continuous records going back 
to the late 1960s; there is thus limited Northern Hemisphere 
mid-latitude information about ozone changes in the 1960s 
and throughout most of the 1970s and virtually no infor- 
mation for other parts of the globe. At the beginning of 
the 1990 high-quality ozone profile measurement at NDACC 
sites started providing reliable long-term information at se- 
lected locations. However, ground-based measurement sta- 
tions are not uniformly distributed around the globe. The 
coverage of the ozone variability at high latitudes and, to 
some extent, in the tropics, were improved by various po- 
lar research projects and with the SHADOZ network. Nev- 
ertheless, there are still large regions with low-measurement 
frequency and sparse coverage. 

More complete information for the global analysis of 
ozone changes only became available with the launch of the 
first satellite instruments. A BUV instrument operated for 
a limited period in the early 1970s with quasi-continuous 
satellite monitoring starting in the late 1970s. The period 
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between the end of the SAGE I measurements and the begin- 
ning of the SAGE II measurements (see Fig. 2) has proven to 
be one of the most challenging obstacles for reliable, highly 
vertically resolved, global ozone trend detection, especially 
with the altitude registration problems reported for SAGE 
I measurements (Wang et al., 1996). SBUV measurements 
are available during this time, however, their coarse verti- 
cal resolution compared to the SAGE measurements severely 
limits their utility in tying the two SAGE records together, 
especially in the lower stratosphere. It is evident that the 
largest ozone decreases in the extra Uopics occurred during 
the 1980s and, therefore, the beginning of this decade is of- 
ten used as the starting point for Uend analyses. However, 
correctly anchoring Uends in this period is complicated and 
represents a challenge for long-term trend analysis (McPeters 
et al., 1994; Randel et al., 1999; Cunnold et al., 2000; Wang 
et al., 2002; Terao and Logan, 2007). 

Another period of sparse global measurements resulted 
from the eruption of Mt. Pinatubo in June 1991. The in- 
creased aerosol loading in the stratosphere resulted in re- 
trieval problems for several satellite instruments (and some 
ground-based instruments) measuring at that time (SAGE II, 
HALOE, SBUV), resulting in recommendations by some in- 
strument science teams that some ozone data for multiple 
months after the Pinatubo eruption be discarded (Yu and She, 
1995; Cunnold et al., 1996; Bhartia et al., 2013). This re- 
duces the available lower stratospheric ozone measurements 
during that period to primarily ground-based instruments 
(ozonesondes). 

One goal of the SI 2 N Initiative is to make the measurement 
community aware of the increased demand for high-quality 
and consistent data sets with well-defined uncertainty esti- 
mates to fulfill requirements for advanced statistical methods 
used to detect small changes in stratospheric ozone variabil- 
ity. It is clear that ODS reductions in response to the Mon- 
treal Protocol has led to a reversal of the negative column 
ozone trends in the extra tropics (e.g. Ziemke and Chandra, 
2012); however, the documentation of the beneficial effect 
on ozone profile changes is still challenging (in particular, 
in an atmosphere where greenhouse gas increases may also 
be playing a role). As stated earlier, this is the first of three 
initiative overview papers and provides a summary of avail- 
able ozone profile measurements. The second overview pa- 
per (Lambert et al., 2014) focuses on comparison and valida- 
tion of the different data sets and the different rettievals for 
the various instruments that are introduced here. A quantita- 
tive understanding of the offsets and drifts among data sets 
and their limitations is essential when merging them for time- 
series analyses. The third overview paper (Harris et al., 2014) 
focuses on describing and characterizing data sets that have 
been derived from a combination of different data sources 
(multiple satellite instruments). Several of these combined 
data sets are compared to each other and a detailed evalua- 
tion of them is given. 


Together, the three overview papers provide the scientific 
community with a source of information about the latest ac- 
tivities of ozone profile measurement groups, with the inten- 
tion of providing needed background information and guid- 
ance for future ozone profile Uend analysis. However, at this 
point it is also important to point out that there remains a sig- 
nificant need for high-quality stratospheric ozone measure- 
ments. Even with the reduction of ODSs in the stratosphere 
following the Montreal Protocol and its Amendments and 
Adjustments an exact attribution of stratospheric ozone re- 
covery is not possible due to changes associated with increas- 
ing greenhouse gas concentrations. Only with highly stable 
and accurate profile ozone measurements will it be possible 
to separate and correctly attribute future ozone changes and 
verify ozone recovery. 
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